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The  energetics,  interfacial  properties,  instabilities,  and  fragmentation  patterns  of  eleetrosprays  made  from 
formamide  salt  solutions  arc  investigated  in  a  mass  speetrometric  vacuum  eleetrospray  experiment  and  using 
molecular  dynamics  (MD)  simulations.  The  eleetrospray  source  is  operated  in  a  Taylor  cone-jet  mode,  with 
the  nanojet  that  forms  being  characterized  by  high  surface-normal  electric  field  strengths  in  the  vicinity  of  1 
V/nm.  Mass-to-eharge  ratios  were  determined  for  both  positive  and  negative  currents  sprayed  from 
Nal— formamide  solutions  with  solute— solvent  mole  ratios  of  1:8.4  and  1:36.9,  and  from  Kl— formamide 
solutions  with  mole  ratios  of  1:41  and  1  83.  The  molecular  dynamics  simulations  were  conducted  on  isolated 
10  nm  Nal— formamide  droplets  at  mole  ratios  of  1:8  and  1:16.  The  droplet  was  subjected  to  a  uniform 
electric  field  with  strengths  ranging  between  0.5  and  1.5  V/nm.  Both  the  experiments  and  simulations 
demonstrate  a  mixed  charge  emission  regime  where  field-induced  desorption  of  solvated  ions  and  charged 
droplets  occurs.  The  macroscopic  parameters,  such  as  average  mass-to-charge  ratio  and  maximum  surface- 
normal  field  strengths  deduced  from  the  simulations  are  found  to  be  in  good  agreement  with  the  experimental 
work  and  consistent  with  electrohydrodynamie  theory  of  cone-jets.  The  observed  mass  speetrometric  Na+ 
and  I"  solvated  ion  distributions  are  consistent  with  a  thermal  evaporation  process,  and  arc  correctly  reproduced 
by  the  simulation  after  incorporation  of  the  different  flight  times  and  unimolccular  ion  dissociation  rates  in 
the  analysis.  Alignment  of  formamide  dipoles  and  field-induced  reorganization  of  the  positive  and  negative 
ionic  charges  in  the  interfaeial  region  arc  both  found  to  contribute  to  the  surface- normal  field  near  the  points 
of  charge  emission.  In  the  simulations  the  majority  of  cluster  ions  arc  found  to  be  emitted  from  the  tip  of  the 
jet  rather  than  from  the  neck  region  next  to  the  Taylor  cone.  This  finding  is  consistent  with  the  experimental 
energy  distributions  of  the  solvated  ions  which  demonstrate  that  indeed  most  ions  are  emitted  closer  to  the 
jet  region,  that  is,  beyond  the  cone-ncek  region  where  ohmic  losses  occur.  This  observation  is  also  consistent 
w  ith  continuum  eleetrohydrodynamic  predictions  of  cluster-ion  evaporation  at  surface  regions  of  high  curvature 
and  therefore  maximum  surface  electric  field  strengths,  which  may  be  the  eonc-neek  region,  the  breakup 
region  of  the  jet  (usually  near  the  tip  of  the  jet),  or  the  emitted  charged  droplets.  In  the  nanoscale  jets  observed 
in  this  study,  the  regions  of  highest  spatial  curvature  are  at  the  ends  of  the  jets  where  nascent  drops  either  are 
forming  or  have  just  detached. 


1.  Introduction 

Understanding  the  fundamental  properties  of  conducting 
liquids  under  the  influence  of  applied  strong  electric  fields  is  a 
topic  of  continuing  and  growing  interest.  The  associated 
phenomenology  governs  a  number  of  venerable  technological 
areas,  most  notably  eleetrospray  ionization1  (see  also  the  Nobel 
lecture  by  John  Fenn  titled  “Eleetrospray  Wings  for  Molecular 
Elephants’'2),  eleetrospray  propulsion,3  eleetrospray  painting, 
inkjet  printing,4  and  electrospinning5  of  materials.  In  all  of  these 
applications,  the  electric  field  causes  deformation  and  charging 
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of  the  liquid  surface,  with  consequent  phenomena  including 
generation  of  jets,  charged-droplet  formation.  Coulomb  fission, 
and  ion-field  evaporation  (also  referred  to  as  field-induced  ion 
desorption).  Most  of  the  above  technological  applications 
involve  breakup  of  the  original  fluid  structure  (often  fluid  jet) 
into  droplets,  with,  possibly,  the  occurrence  of  subsequent 
disintegration  processes.  Indeed,  we  note  here  that  the  use  of 
electric  fields  for  controlling  fluids  dates  back  to  (at  least)  the 
eighteenth  century,  with  demonstrations  of  the  influence  of 
strong  electric  fields  on  the  deformation  and  disintegration 
patterns  of  water  jets  appearing  already  in  1749,6  and  the 
employment  of  electric  fields  for  the  deflection  of  liquid  jets 
described  in  1859.7  The  breakup  of  a  charged  liquid  drop,  also 
called  Coulomb  fission,  was  analyzed  in  1882  by  Rayleigh,  who 
showed  that  such  instability  will  occur  when  the  forces  due  to 
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electrostatic  (Coulombic)  repulsion  of  surface  charges  (expressed 
as  the  electrostatic  pressure  for  a  charged  sphere,  Pe  =  <J2/2f:<), 
where  o  is  the  surface  charge  density  and  is  the  permittivity 
of  vacuum)  exceed  the  forces  associated  with  the  surface  tension 
of  the  liquid  (expressed  as  the  pressure  associated  with  the 
curvature  and  surface  tension  y  of  a  liquid  spherical  droplet  of 
radius  R,  Pr  =  2  y/R).  The  condition  where  the  two  forces  are 
equal  is  referred  to  as  the  Rayleigh  instability  limit.8  Rayleigh 
showed  that  the  charge  of  a  droplet  at  this  limit,  with  a  diameter 
1 1 ,  is  given  by 

qR  =  (&jT£0yd})U2  (1) 

The  electric  held  strength  at  the  droplet  surface  in  the  absence 
of  an  external  held  is  given  by 

£  =  -2-;  (2) 

m0d- 

and  the  critical  held,  £R,  is  obtained  from 


The  first  observation  of  the  “Rayleigh  instability'’  was 
reported  by  Zeleney  about  thirty  years  after  the  work  of 
Rayleigh.9-10  In  these  pioneering  experiments,  Zeleney  passed 
a  stream  of  conducting  volatile  liquid  through  a  small  bore,  thin- 
walled  tube,  maintained  at  a  high  potential  relative  to  an 
opposing  counter  electrode.  The  resulting  intense  electric  held 
at  the  tip  of  the  tube  dispersed  the  emanating  liquid  into  a  hne 
spray  of  charged  droplets.  Zeleney  was  able  to  observe  the 
breakup  of  the  charged  droplets  as  the  solvent  evaporated  (with 
the  remaining  droplet  thus  reaching  the  stability  limit,  as 
predicted  by  Rayleigh).  Seventy  years  later,  Gomez  and  Tang 
published  some  striking  photographic  records  of  electrospayed 
heptane  drops  undergoing  a  Coulomb  explosion.11  12 

We  note  here  that  although  the  above  prediction  for  q R  has 
been  reported  to  agree  approximately  with  most  measurements,14 
more  rehned  observations  indicate  that  the  instability  often 
occurs  at  about  80%  of  r/R.14~16  The  generation  of  a  fine  spray 
via  the  application  of  an  electric  field,  as  in  Zeleney's  experi¬ 
ments,  is  a  consequence  of  the  field-induced  ‘focusing'  of  fluids 
into  pointed  structures,  referred  to  genencally  as  ‘Taylor  cones”; 
Sir  Geoffrey  Ingram  Taylor  was  primarily  interested  in  the 
behavior  of  water  droplets  in  strong  electric  fields,  such  as  in 
thunderstorms.17  This  cone  geometry  corresponds  to  the  hy¬ 
drostatic  structure  in  which  a  balance  between  the  capillary  and 
electrostatic  stresses  is  exactly  satisfied  on  the  surface  of  an 
equipotential  cone.  Development  of  the  Laplace  equation  in 
polar  coordinates  leads  to  the  condition  P t/2(cos  0)  =  0  on  the 
L  egendre  function  of  order  l/2,  with  the  solution  0  =  49,3°  (often 
referred  to  as  the  ‘Taylor  angle”)  for  the  cone  half  angle  for 
which  the  above  stress  balance  is  obtained  on  the  conical 
equipotential  surface.  For  a  comprehensive  review  of  the  current 
literature  on  the  fluid  dynamics  of  Taylor  cones,  see  a  recent 
account  by  Fernandez  de  la  Mora.18 

Although  Taylor’s  treatment  applied  to  a  nonconducting  liquid 
or  neutral  molecules  with  dipole  moments  (for  which  there  is 
no  electric  current  and  no  liquid  flow),  his  name  continues  to 
be  used  also  in  applications  where  the  liquid  conducts  (due  to 
the  presence  of  charged  species,  i.e.,  cations  and  anions).  In 
the  latter  case,  a  thin  jet  may  emerge  from  the  tip  of  the  cone, 
and  it  will  break  up  into  a  series  of  droplets  because  of  another 
instability  described  first  by  Rayleigh  for  uncharged  liquids19-20 
called  a  “varicose  instability”;  this  instability  has  been  observed 
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by  Plateau  in  1872. 21  Underlying  this  instability  are  surface 
tension  and  capillary  effects  that  are  expressed  as  waves  that 
grow  in  amplitude  until  a  droplet  is  pinched  off.  The  accepted 
mechanism  of  jet  breakup  involves  flow  from  regions  of  the 
liquid  column  (cylindrical  jet  or  bridge)  with  smaller  radii,  that 
is,  regions  of  higher  curvature  (i.e.,  necks  where  the  Laplace 
pressure  is  larger)  to  neighboring  crest  regions  with  lower 
curvature  (and  hence  lower  pressure).  Because  the  resulting 
curvature  in  the  necks  becomes  larger  and  larger  as  the  liquid 
flows  away  from  this  region,  the  flow  of  liquid  from  the 
shrinking  neck  to  the  bulging  crest  accelerates,  leading  to  a 
singularity  in  finite  time  (namely,  the  flow  velocity  diverges), 
which  is  accompanied  by  pinch-off  of  a  droplet  (see  ref  69  and 
references  cited  therein). 

Two  scenarios  for  the  subsequent  evolution  of  droplets 
generated  via  jet  breakup  have  been  proposed.  The  first  is  the 
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one  proposed  by  Dole  and  co- workers.22  Dole’s  charged  residue 
model  (CRM)  has  been  developed  to  explain  experiments  on 
electrically  driven  jets  of  a  solution  of  polystyrene  molecules, 
generated  by  a  method  (similar  to  that  used  earlier  by  Zeleney) 
which  Dole  described  in  his  biography23  as  “1  got  this  idea  from 
learning  about  the  electrospaying  of  paint  onto  automobile 
bodies  while  working  as  a  consultant  for  a  paint  company  in 
Chicago”.  In  the  CRM,  large  charged  droplets  of  volatile  liquids 
will  evaporatively  diminish  in  size  (that  is,  reduce  their 
diameters)  while  maintaining  the  charge,  ultimately  reaching 
the  Rayleigh  limit  and  undergoing  Coulomh  fission;  this  process 
is  assumed  to  repeat  itself  until  a  charged  solute  residue  is 
formed. 

An  alternative  model  (called  the  ion  evaporation  model,  or 
1EM)  for  the  evolution  of  a  large  charged  droplet  subsequent 
to  its  generation  via  jet  breakup,  has  been  advanced  by  Iribame 
and  Thomson,24  two  meteorologists  who  were  interested  in  the 
possibility  that  charged  droplets  of  seawater  might  be  a  source 
of  some  of  the  ions  found  in  the  atmosphere.  These  researchers 
explored  the  conditions  at  which  the  electric  held  strength  at 
the  surface  of  the  droplet  exceeds  a  threshold  value,  £j,  at  w  hich 
ion  held  evaporation  occurs.  This  condition  can  be  met  provided 
£i  <  ER(d)  (otherwise  Coulomb  fission  occurs),  and  the  solute 
residue  diameter  is  less  than  the  critical  diameter  at  which  E| 
=  Er (d)\  i.e.,  in  the  case  of  salt  solutions  the  concentration  is 
low.  The  threshold  value  E\  derives  from  the  energy  to  overcome 
the  barrier  associated  with  the  combination  of  the  solvation 
energy  of  the  evaporated  ion,  and  the  electrostatic  interaction 
between  the  ion  and  the  droplet,  lnbame  and  Thomson24 
estimated  that  this  barrier  is  approximately  4—5  A  from  the 
droplet  surface. 

On  the  basis  of  a  kinetic  model  and  charge  mobility 
experiments,  Thomson  and  lnbame25  concluded  that  the  critical 
held,  Ei  was  approximately  1  V/nm,  and  that  field  evaporation 
would  be  expected  for  droplets  with  diameters  <  10  nm. 
Loscertales  and  Fernandez  de  la  Mora26  and  Katta  and  co¬ 
workers27  used  charge  mobility  experiments  and  refined  models 
to  determine  electric  fields  between  0.7  and  1.0  V/nm  for  the 
evaporation  of  different  singly  charged  ions.  Fenn  et  al.28  also 
arrived  at  similar  fields  by  applying  a  simple  electrostatic  model 
to  the  interpretation  of  their  experiments  on  dilute  polyethylene 
glycol  solutions,  assuming  that  the  charges  of  small  droplets 
are  equally  spread  on  the  droplet  surface,  and  that  the  charge 
distances  on  the  droplet  cannot  be  less  than  that  found  on  the 
multiply  charged  large  field-evaporated  ions.  These  authors 
determined  fields  ranging  from  1,46  to  2.66  V/nm  with 
increasing  field  strength  required  for  smaller,  less  charged 
oligomer  ions. 

The  investigation  of  the  detailed  field-evaporation  mecha¬ 
nisms  from  charged  droplets  is  most  challenging  because  a 
charge-shedding  isolated  droplet  does  not  ever  attain  a  steady 
state.  The  field-induced  evaporation  phenomenon,  however,  is 
not  limited  to  droplets.  Indeed,  the  conditions  for  field  evapora¬ 
tion  can  also  be  met  at  the  surface  of  a  flowing  electrified  jet. 
The  experimental  detection  and  analysis  of  nascent  ion  currents 
from  jets  requires  spray  operation  under  vacuum  conditions 
where  conventional  ion  optics  can  be  used  for  the  analysis  of 
charged  particles  prior  to  collisions  with  a  background  gas. 
Because  the  use  of  a  volatile  solvent  such  as  water  or  methanol 
in  a  vacuum  would  lead  to  rapid  freezing  of  the  jet,  the  study 
of  field-evaporation  from  a  jet  is  limited  to  liquids  with  low 
vapor  pressures. 

In  Figure  1  we  display  schematically  different  configurations 
of  electrosprays  involving  a  capillary  and  an  extraction  electrode. 


pure  ion  mixed  ion-  pure  droplet 
droplet 


Figure  1.  Schematic  cross-sectional  representation  of  several  con¬ 
figurations  and  ejection  modes  of  a  capillary  elecirospray.  The  extraction 
voliage,  V,  is  applied  between  an  extraction  electrode  (shaped  as  u  ring) 
and  the  capillary  containing  the  conducting  liquid  (a  tormamide  sail 
solution  in  the  preseni  study).  Vc  is  the  critical  extraction  voliage  ai 
which  a  Taylor  cone  is  formed.  The  lop  row  illustrates  ihe  liquid 
response  as  a  function  of  the  extraction  voliage:  from  field-free 
conditions,  to  a  stressed  jet  operation,  V  >  Vc.  The  bottom  row  illustrates 
the  different  mixed  modes  of  operation,  when  the  jel  is  operated  in  a 
single  cone-jet  configuration  (V  ^  Vc).  The  darkly  shaded  area  outside 
the  orifice  of  ihe  capillary  represents  ion  emission;  this  is  shown  for 
ihe  case  of  pure  ion  emission  (leftmost  configuration  in  ihe  bottom 
row),  as  well  as  for  a  mixed  ion-droplet  emission  mode  (middle 
configuration  in  the  boiiom  row). 

The  study  of  held  evaporation  on  jets  is  best  accomplished  at 
well  characterized  conditions  provided  by  extraction  voltages 
close  to  the  critical  voltage,  V,  at  which  the  Taylor  cone  that 
triggers  electrospray  is  formed.  Depending  on  the  choice  of 
liquid,  operation  in  this  cone-jet  configuration  can  lead  to  three 
different  modes  of  charge  emission:  a  pure  ion  emission  mode, 
a  mixed  ion-droplet  emission  mode,  and  a  pure  charged  droplet 
emission  mode.  Liquid  metals,  having  very  high  surface 
tensions,  and  consequently  high  resistance  to  Coulomh  fission, 
arc  known  to  produce  pure,  intense  metal  ion  beams  in  an 
electrospray  configuration  with  little  or  no  contribution  from 
charged  droplets.29  Jet  formation  does  not  occur  under  typical 
liquid  metal  electrospray  conditions,  and  the  Taylor  cone  is 
stabilized  through  pure  ion  emission  at  the  tip  of  the  Taylor 
cone  where  the  surface  normal  fields  are  highest.  Solutions 
involving  low-vapor  pressure  solvents,  which  are  the  subject 
of  this  work,  have  significantly  more  flexible  electrospray 
properties.  Such  systems  generally  operate  in  a  mixed  ion- 
droplet  emission  regime  with  complex  mass  spectra  of  the 
emitted  ions  that  accompany  charged  droplets.  At  normal  cone- 
jet  conditions,  with  the  exception  of  work  discussed  further 
below,  the  ion  component  represents  a  very  minor  fraction  of 
the  total  emission  current.  Pure  charged  droplet  emission  has 
only  recently  been  demonstrated30  by  the  suppression  of  ion 
field  evaporation  through  the  choice  of  solutes  with  very  high 
solvation  energies. 

Most  of  the  early  work  on  electrosprays  of  organic  solutions 
in  vacuum  was  conducted  with  glycerol.  Huberman31  conducted 
an  energy  analysis  of  positive  ions  sprayed  from  a  Nal -glycerol 
solution,  and  detected  a  feature  that  was  attributed  to  metastable 
solvated  ions  losing  a  solvent  molecule.  Stimpson  et  al.32 
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Figure  2.  Schematic  representation  of  an  electrospray  in  the  cone-jet 
regime  (see  Figure  1  for  V  ~  Vc).  E  represents  the  surface  normal 
electric  field  and  Q  is  the  volume  flow  rate.  Note  the  varicose  wave  on 
the  surf  ace  of  the  jet.  The  electric  held  at  the  neck  (E)  and  the  external 
field  (£appUcd)  are  denoted.  In  the  mixed  ion-droplet  regime  (see  the 
middle  configuration  in  the  bottom  row  of  Figure  1 )  the  ejected  spray 
consists  of  a  mixture  of  droplets  and  ions  (the  latter  not  shown  here). 


observed  similar  behavior  involving  metastable  species  in  a  mass 
spectrometric  analysis  of  both  positive  and  negative  ions  emitted 
from  glycerol  solutions  of  a  number  of  salts.  A  large  number 
of  solvated  ions  as  well  as  protonated  glycerol  ions  were 
detected.  In  the  experiments  by  Huberman  and  Stimpson  et  al., 
very  high  extraction  voltages,  exceeding  8  kV,  were  used.  The 
high  resulting  fields  are  known  to  generate  a  stressed  jet  regime, 
leading  to  multiple  jets  on  the  rim  of  the  capillary  orifice  (see 
top  right  in  Figure  1  ).33  In  recent  years,  notable  advances  have 
been  made  in  understanding  electrosprays  operating  in  the 
significantly  more  characterized  single  cone-jet  mode,34-43 
depicted  in  greater  detail  in  Figure  2.  As  aforementioned, 
Taylor17  showed  that  the  surface  tension  and  electrostatic  forces 
for  a  nonconducting  liquid  can  balance  in  a  way  that  leads  to 
formation  of  a  liquid  cone.  Above  a  critical  field,  the  cone 
becomes  unstable  and  a  jet  is  formed;  conditions  for  formation 
of  the  cone-jet  mode  have  been  discussed33  early  on  and  most 
recently  Ganan-Calvo44  described  the  cone -jet  mode  as  a 
combination  of  a  slender  jet  and  a  static  Taylor  cone  near  the 
capillary  (nozzle)  orifice.  Fernandez  de  la  Mora  derived  a  more 
general  theory  for  single  cone-jet  sprays  of  conducting  liquids 
predicting  a  broader  range  of  cone  apex  angles  w  hich  was  also 
confirmed  experimentally.34  The  extraction  voltages  to  operate 
in  such  a  Taylor  cone-jet  mode  are  much  lower  (on  the  order 
of  1  —2  kV)  than  those  used  in  the  earlier  experiments,  leading 
to  a  significantly  gentler  break-up  of  the  jet.  In  ground-breaking 
work,  Fernandez  de  la  Mora  and  Loscertales35  derived  both 
theoretically  and  experimentally  an  expression  that  gives  the 
proportionality  between  the  parameters  of  a  conducting  liquid 
and  the  current  emitted  by  the  Taylor  cone 

1~(YKQ/e)'p-  (4) 

where  K  and  e  are  the  liquid  conductivity  and  dielectric  constant, 
respectively,  and  Q  is  the  volume  flow  rate.  It  is  noteworthy 
that  no  electrostatic  parameters  appear  in  this  expression.  The 
role  of  the  external  field  is  the  generation  of  the  Taylor  cone 
and  the  associated  surface-charge  buildup.  Equation  4  can  be 
rearranged  and  rewritten  in  terms  of  the  specific  charge  of 
the  emitted  droplets,  qfm , 

4_i  fiE.  (5) 

m  pV  rQ  <5) 

where  p  is  the  liquid  density.  Because  droplets  formed  in  jet 
fission  are  expected  to  be  charged  near  the  Rayleigh  limit,  the 
specific  charge  is  also  a  measure  of  the  diameter,  Thus,  eq  5 
states  that  the  droplet  size  emitted  by  a  Taylor  cone  can  be 


controlled  hy  parameters  such  as  the  liquid  conductivity  and 
the  volume  flow  rate. 

Inspection  of  the  cone-jet  structure  in  Figure  2  allows  the 
identification  of  three  regions:  the  cone,  the  jet  and  the 
intermediate  neck  region.  The  last  represents  the  transition  from 
a  flow  region  in  the  Taylor  cone,  in  which  the  electrostatic 
nonequilibnum  associated  with  the  narrowing  cone  diameter 
and  increasing  flow  speed  leads  to  ohmic  currents  within  the 
liquid,  to  a  convective  flow  region  in  the  jet  that  is  in 
electrostatic  equilibrium,  and  where  the  charge  density  is 
lowered  as  the  jet  speed  increases.  Thus,  the  region  of  highest 
surface  charge  and  normal  electric  field  is  in  the  intermediate 
neck  region.  Gamero-Castano  and  Femdndcz.  dc  la  Mora41 
derived  expressions  for  the  maximum  normal  electric  field  and 
radius  of  curvature  or  jet  radius  in  this  region: 


M)  ^ 

Km,x  =  [e0Q/K]w 


(6) 

(7) 


The  above  authors  were  also  the  first  to  observe  field  evaporation 
from  a  jet  using  Nal— formamide  solutions.41  In  these  experi¬ 
ments,  volume  flow  rates  of  less  than  100  pl/s  and  conductivities 
exceeding  1  S/m  were  required,  resulting  in  electric  fields  near 
1  V/nm,  which  are  consistent  with  earlier  predictions  for  the 
electric  field  values  that  correspond  to  the  generation  of  field- 
induced  ion  desorption  events.  Conditions  were  achieved  at 
which  the  ion  current  exceeded  the  charged  droplet  current.  At 
these  more  balanced  mixed  ion/droplet  emission  conditions,  the 
curvature  of  the  transition  region  and  the  diameter  of  the  jet 
are  predicted  to  be  in  the  vicinity  of  10  nrn.  This  is  again  in 
full  analogy  to  the  electrospray  ionization  mechanism  postulated 
for  droplets. 

Subsequently,  ion  field  evaporation  has  been  commonly 
observed  with  a  third  class  of  low-vapor  pressure  liquids, 
namely,  ionic  liquids45-47  (the  first  two  being  liquid  metals  and 
the  organic  liquids  noted  above).  Very  low  flow  rates  can  be 
obtained  by  replacing  the  capillary  with  a  needle  etched  to 
optimize  wetting  properties.48  Pure  ion  emission,  signifying  a 
vanishing  jet,  has  now  been  observed  for  the  ionic  liquids 
[Emim)[BF4|47  and  [Emim|[ Iin]  (Emim,  l-ethyl-3-methylim- 
idazolium;  Im,  b/.v(trifluoromethylsulfonyl)imide),49  whereby  the 
former  has  been  successfully  sprayed  using  a  microfabricated 
needle  array.50  In  the  present  paper  we  focus  on  formamide  salt 
solutions  because  of  the  rich  variety  of  field  evaporated  tons 
that  are  found  for  these  systems,  and  because  they  are  more 
computationally  tractable,  thus  allowing  systematic  comparisons 
between  the  experimental  measurements  and  theoretical  predic¬ 
tions.  Moreover,  as  demonstrated  in  several  electrospray  ioniza¬ 
tion  studies,  ion  field  evaporation  at  fields  near  the  onset  of  ion 
desorption  may  be  viewed  as  a  soft-ionization  approach,  where 
the  properties  of  ions  produced  in  this  “gentle*’  desorption 
process  reflect  the  microscopic  properties  of  the  parent  solution 
and  its  interface,  as  was  assumed  in  the  original  analysis  hy 
Fenn  and  co-workers.28 

Recently,  Chiu  et  al,4"  conducted  mass  spectrometric  analysis 
of  positive  ions  emitted  from  a  concentrated  (28.4  wt  %) 
Nal— formamide  solution  operated  in  a  Taylor  cone-jet,  ion  field 
evaporation  mode.  The  measurements  identified  primarily 
solvated  ions,  Na+(HCONFh),„  with  a  preferred  solvation 
number  of  n  =  4  (the  number  of  solvating  formamide  molecules 
per  solvated  Na+),  which  is  the  largest  solvated  ion  predicted 
to  be  thermally  stable  with  respect  to  the  experimental  ion  flight 
times  at  room  temperature.  Thus  the  measurements  appear  to 
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confirm  the  notion  of  a  thermal  evaporation  process.  Interest¬ 
ingly,  however,  the  solvated  ion  energy  distributions  were  found 
to  be  near  or  below  the  energies  of  the  charged  droplets, 
contrasting  the  electrohydrodynamic  prediction  that  ion  held 
evaporation  would  occur  from  the  neck  region  of  the  jet  which 
is  expected  to  have  a  higher  potential  than  the  tip  of  the  jet 
because  of  ohmic  losses  in  the  cone -jet  transition  region.  Similar 
observations  were  made  by  Gamero-Castano43  for  the  same 
solution.  However,  with  measurements  on  the  ionic  liquid 
[Emim][Im]  in  a  mixed  ion/droplet  regime,  Chiu  et  al.51-52  could 
demonstrate  clearly  that  field-evaporated  ions  are  emitted  at  the 
cone-jet  neck,  as  evidenced  by  their  energy  and  angular 
distribution,  consistent  with  theoretical  arguments.41 

Recent  experimental  and  theoretical  advances53  55  have  shed 
new  light  on  the  interfacial  properties  of  aqueous  salt  solutions 
under  conditions  free  of  external  electric  fields.  Counter  to  the 
long-held  view  that  atomic  ions  are  repelled  by  the  air/water 
interface  into  the  bulk,  MD  simulations  of  alkali-halide  solutions 
have  demonstrated  that  more  polarizable  anions  reside  prefer¬ 
entially  at  the  liquid-to- vacuum  interface.53-  59  These  results  have 
been  confirmed  by  X-ray  photoelectron  spectroscopy  experi 
ments,60  as  well  as  sum-frequency56-57-61*62  and  second-harmonic 
generation  experiments63-67  of  aqueous  halide  solutions  at 
ambient  pressure.  Though  these  studies  examine  the  microscopic 
hehavior  of  static  solution  interfaces,  MD  simulations68  have 
heen  conducted  on  an  acidic  positively  charged  water  nano- 
droplet,  with  the  amount  of  excess  charge  near  the  Rayleigh 
limit.  In  these  simulations,  ion  held  evaporation  of  hydronium 
ions  was  observed;  we  note  here  that  in  these  simulations 
negative  ions  were  not  included,  and  this  may  influence  the 
properties  of  the  simulated  system,  particularly  at  reduced 
dimensions. 

We  conclude  our  brief  review  of  recent  work  with  MD 
simulations,  performed  by  Moseler  and  Landman,69  of  liquid 
flows  at  nanoscale  dimensions,  and  the  formulation  of  a 
stochastic  hydrodynamic  treatment.  In  these  studies,  simulations 
were  performed  of  liquid  propane  flow  ing  under  applied  back¬ 
pressure  through  a  cylindrical  tube,  and  emanating  through  a 
nanoscale  nozzle  into  a  vacuum;  recently,  Kang  and  Landman 
have  extended  these  simulations  to  injection  of  liquid  jets  into 
a  surrounding  fluid  70  The  properties  of  the  resulting  jets  —  in 
particular  the  characteristics  of  the  breakup  of  the  liquid  jet  into 
droplets  —  were  found  to  differ  from  continuum  hydrodynamics 
predictions.  It  was  discovered  that  a  newly  derived  hydrody¬ 
namic  formulation  that  includes  intrinsic  size-dependent  fluctua¬ 
tions,  provides  a  faithfully  accurate  descnption,  in  agreement 
with  the  atomistic  simulations.  This  work  lays  the  groundwork 
for  an  investigation  of  electrospray  jets,  where  the  cone-jet 
transition  region  can  be  regarded  as  a  virtual  nozzle  of  nanoscale 
dimensions. 

In  this  paper  we  present  a  systematic  study  of  ion  field 
evaporation  from  salt  solutions  using  a  combination  of  mass 
spectrometry  and  MD  simulations.  To  gain  insights  into  the 
microscopic  processes  of  nanojet  formation  and  charge  emission 
from  dielectric  liquids  under  the  influence  of  strong  electric 
fields,  we  chose  to  investigate  theoretically  the  behavior  and 
response  of  liquid  drops  (made  of  pure  formamide  or  Nal  — 
formamide  solutions)  because  studies  on  such  systems  free  us 
from  complicating  factors,  including  considerations  of  effects 
due  to  flow  rates,  contact  angles,  and  the  geometry  of  the 
capillary  tube  or  needle,  used  to  generate  the  slender  liquid 
flowing  system;  for  a  discussion  of  the  analogy  between  cone- 
jets  on  an  exploding  drop  and  those  forming  on  a  meniscus 
supported  by  a  capillary  tube,  see  the  article  by  Fernandez  de 


la  Mora16  entitled  “On  the  Outcome  of  the  Coulombic  Fission 
of  a  Charged  Isolated  Drop”.  The  MD  simulations  were 
conducted  on  a  10  nm  diameter,  uncharged  (equal  numbers  of 
positive  Na+  and  negative  1  charges)  Nal— formamide  droplet 
subjected  to  external  electric  fields  with  strengths  on  the  order 
of  1  V/nm.  The  size  of  the  drop  was  selected  such  that  the 
corresponding  curvature  of  the  droplet  (1  /R)  would  be  compa¬ 
rable  to  the  estimated  minimum  curvature  expected  to  be 
encountered  on  a  liquid  jet  (see  eq  7),  and  thus  the  physics  of 
jet  formation,  ion  field  evaporation,  and  charge  emission  from 
such  drops  can  be  taken  as  a  faithful  representation  of  the 
processes  occurring  in  electric  field-driven  cone-jet  flow  systems. 
Indeed,  experimental  results  of  field-induced  dual-jet  formation 
on  droplets,  either  flowing  through  a  region  of  high  electric  field 
or  levitated,  have  been  reported;71-74  these  studies  relied  on 
visual/optical  examination  and  were  conducted  on  droplets  with 
diameters  exceeding  50  /tm  exhibiting  formation  of  jets.  Under 
these  conditions,  ion  field  evaporation  from  the  main  drop  is 
not  expected  to  occur.  We  remark  here  that  Grimm  and 
Beauchamp  conducted  mass  spectrometric  measurements  of  field 
evaporated  ions  from  droplets  of  a  methanol  solution  that  are 
ejected  from  larger  (170  /<m  diameter)  drops.74  However,  the 
high  external  fields  that  would  be  required  to  induce  direct  field 
desorption  from  10  nm  diameter  droplets  (which,  as  noted  above, 
is  the  size  range  of  interest,  here)  were  outside  the  range  of 
these  earlier  experiments. 

The  mass  spectrometric  measurements  presented  in  this  paper 
were  conducted  on  ions  that  are  field  evaporated  from  a  cone- 
jet  of  solutions  (Nal  and  KI  dissolved  in  formamide)  with  similar 
concentrations  to  those  examined  in  the  simulations.  Quantum 
chemical  calculations  were  carried  out  to  determine  the  energy- 
optimal  structures  and  thermochemistry  of  solvated  ion  clusters 
with  sizes  equal  to  those  observed  experimentally.  The  ther¬ 
mochemistry  and  associated  normal-mode  frequencies  are  used 
to  assess  the  thermal  stabilities  and  lifetimes  of  the  ions.  The 
solvated  ion  cluster  lifetimes  are  of  particular  importance  for 
our  analysis  because  they  are  used  in  extrapolating  and 
comparing  the  simulated  results  concerning  mass  distributions 
of  emitted  solvated  ions  (that  are  simulated  for  relatively  short 
flight  times)  to  the  experimental,  mass  spectrometric  data, 
obtained  for  much  longer  flight  times  than  the  simulated  ones. 
Both  positive  and  negative  solvated  ion  distributions,  obtained 
experimentally  from  respective  positive  and  negative  extraction 
fields,  and  from  calculated  evaporation  statistics  at  opposite  ends 
of  the  simulated  droplet,  are  presented  and  analyzed. 

2.  Experimental  Section 

The  mass  spectrometric  experiment  has  been  described  in 
detail  previously,47  and  only  the  features  relevant  for  this  article 
will  be  discussed.  A  schematic  of  the  apparatus  is  shown  in 
Figure  3.  The  electrospray  source  consists  of  a  tapered  5  /mi 
inner  diameter  fused  silica  capillary  tip  (New  Objective,  Inc.) 
that  is  fed  by  a  50  //m  inner  diameter  transport  capillary  via  a 
capillary  union  (Upchurch  Scientific).  The  solution  flow  rate  is 
controlled  by  pressurizing  the  liquid  reservoir  supplying  the 
transport  capillary  with  dry  nitrogen.  The  flow  rate  is  determined 
from  the  velocity  of  bubbles  introduced  into  the  liquid  near  the 
transport  capillary  entrance.  The  capillary  tip  is  positioned  within 
1  mm  of  an  extractor  electrode  with  an  orifice  of  6  mm  diameter. 
The  tip  and  extraction  optics  are  all  in  a  vacuum  (10“7  Torr). 
Typical  extractor  voltages  used  for  the  28.4  wt  %  Nal -formamide 
solution  were  2  kV.  The  capillary  bias  voltage  was  applied  to 
the  capillary  through  a  metal  layer  that  coated  both  the  inner 
and  outer  surface  of  the  capillary  at  the  entrance  end.  Typical 
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Figure  3.  Schematic  diagram  of  the  mass  spectrometnc  experiment 
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bias  potentials  were  550  V,  whereas  the  extractor  was  biased 
at  —1450  V  for  positive  polarity,  and  the  reverse  for  negative 
polarity,  i.e.,  —550  and  +1450  V.  The  capillary  bias  potential 
was  chosen  to  minimize  the  ion  energy  in  the  mass  filter,  which 
is  biased  at  ground  potential,  while  maintaining  reasonable  ion 
transmission  through  the  apparatus. 

The  charged  particle  beam  extracted  from  the  capillary  tip 
can  be  focused  either  on  a  near-field  target  consisting  of  an 
interchangeable  Faraday  cup,  or  a  quartz  crystal  microbalance 
(QCM;  XTM/2,  Inficon),  or  toward  the  entrance  of  a  quadrupole 
mass  filter.  The  near-field  target  is  used  for  macroscopic 
characterization  of  the  charged-particle  beam  through  total 
current  and  mass  flow  measurements.  The  ratio  between  mass 
(low  and  current  provides  the  average  mlq  value,  which,  in 
conjunction  with  ion  mass  distribution  measurements,  provides 
a  measure  of  the  fraction  of  charged  droplets  in  the  beam  The 
near- field  target  is  translated  away  from  the  beam  axis  for  mass 
spectral  measurements.  Given  the  possibility  of  broad  ion  energy 
distributions,  the  quadrupole  filter  is  the  mass  spectrometer  of 
choice  because  the  quadrupole  mass  setting  is  independent  of 
energy.  Fiowever,  there  is  a  significant  penalty  in  mass  resolution 
that  declines  with  ion  energy.  Following  transmission  through 
the  quadrupole,  the  ions  pass  through  a  retarding  potential 
analyzer  prior  to  being  detected  using  a  channel  electron 
multiplier.  The  retarding  potential  analyzer  is  used  to  measure 
the  energy  distribution  of  mass  selected  ions. 

Mass  calibration,  the  mass  transmission  function  and  detector 
response  are  obtained  from  mass  spectra  of  a  standard  calibrant 
compound,  FC-43,  with  a  known  fragmentation  pattern  at  an 
electron  energy  of  70  eV.  The  spectra  were  obtained  after 
replacing  the  electrospray  source  with  an  electron-impact 
ionization  source.  The  intensities  of  fragmented  ions  of  FC-43 
were  measured  for  all  investigated  quadrupole  transmission 
energies,  from  which  the  mass-dependent  discrimination  cor¬ 
rection  factors  were  determined. 

The  concentrated  28.4  wt  °fc  and  1/4  volumetrically  diluted 
Nal  —  formamide  solutions  used  in  this  study  have  salt-to- 
formamide  mole  ratios  of  1:8.4  and  1:36.9,  respectively. 
Similarly,  the  salt/formamide  mole  ratios  determined  for 
Kl— formamide  solutions  in  this  study  are  1:40.9  and  1:83.0, 
respectively.  The  measured  conductivities  for  the  concentrated 
and  dilute  Nal— formamide  solutions  are  1.64  and  1.26  S/m, 
respectively;  conductivities  of  1.32  and  0.77  S/m  are  found  for 
the  concentrated  and  dilute  KI— formamide  solution. 

3.  Computational  Methods 

3.1.  Molecular  Dynamics  Simulations.  To  enable  large- 
scale  atomistic  simulations  of  the  thermodynamic  and  electro¬ 
hydrodynamic  properties  of  complex  fluids  (with  large  partial 
and  ionic  charges)  under  the  influence  of  strong  electric  fields, 
a  computer  code  with  parallel  architecture  was  developed.  To 


model  these  sy  stems  faithfully,  it  is  important  to  treat  the  long- 
range  Coulombic  interactions  accurately  and  efficiently.  In  our 
initial  studies  of  pure  formamide  droplets,  long-range  Coulombic 
interatomic  potentials  were  treated  with  finite  interaction  cut¬ 
offs  ranging  from  1  to  4  nm.  We  found  that  the  employment  of 
finite  range  cutoffs  inhibited  the  expected  droplet  elongation 
as  predicted  by  continuum  theory74'  of  dielectric  drops.  Conse¬ 
quently,  we  concluded  that  Coulombic  interactions  must  be 
treated  accurately  (namely  with  no  cutoff)  to  bring  the  simulation 
results  into  (at  least)  qualitative  agreement  with  existing 
continuum  theory.  Obviously,  the  addition  of  ionic  species 
further  increases  the  necessity  of  proper  implementation  of  long- 
range  Coulomb  interactions.  Indeed,  the  more  accurate  treatment 
of  dipole— dipole  and  dipole-field  cooperative  interactions  allows 
droplets  to  respond  correctly  to  external  fields,  including  field- 
induced  polarization,  elongation,  and  instabilities  that  result  in 
emission  of  charged  cluster  ions  and  electrified  jets, 

In  our  current  simulations  the  long-range  forces  are  modeled 
through  the  use  of  the  fast  multipole  method  (FMM)76  specif¬ 
ically  formulated  and  programmed  for  a  parallel  architecture 
and  adapted  to  our  present  system,  consisting  of  an  isolated 
droplet.  The  FMM  algorithm  enables  the  modeling  of  large 
systems  with  long-range  interactions  by  using  a  “tree- structure” 
subdivision  of  cells  that  allows  a  charge  to  interact  with  a  small 
set  of  multipoles  representing  entire,  increasingly  larger,  volumes 
of  remote  space,  rather  than  with  all  of  the  individual  charges 
contained  within  the  cells.  Important  features  of  the  program 
developed  for  this  study  are 

(1)  The  ability  to  replicate  basic  root-cells  and  “stack”  them 
along  a  particular  axis  to  create  a  long  narrow  tree-structure. 
This  is  particularly  useful  for  efficient  simulations  of  systems 
such  as  highly  elongated  droplets  and  long  liquid  jets. 

(2)  A  large  degree  of  flexibility  in  the  choice  of  boundary 
conditions  (bes)  The  computer  program  allows  the  implementa¬ 
tion  of  periodic  boundary  conditions  (pbes)  in  one,  two  or  three 
spatial  directions,  or  the  execution  of  “cluster  simulations”  where 
pbes  are  not  imposed  altogether.  In  simulations  using  pbes  in 
less  than  three  directions,  the  nonperiodic  directions  can  be 
reflecting  (closed  systems)  or  absorbing  (open  system).  Absorb¬ 
ing  bes  (any  atoms  that  cross  the  system  boundary  are  removed 
from  the  simulation)  tire  necessary  in  simulations  of  fragmenting 
droplets  or  liquid  jets  where  there  are  large  fluxes  of  atoms 
crossing  the  finite  boundaries  (keeping  these  atoms  in  the  root¬ 
cell  would  cause  artificial  effects). 

(3)  A  modified  FMM  algorithm  so  that,  when  using  pbes, 
the  tree- structure  of  the  root-cell,  along  with  all  its  computed 
multipoles,  can  be  replicated  a  specified  number  of  times  as 
neighboring  image  cells.  When  this  is  done,  atoms  in  the  central 
computational  cell  see  the  correct  multipole  structure  as  they 
interact  with  neighboring  periodic  image  cells  out  to  any  desired 
distance. 
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(4)  The  ability  to  perform  constant  pressure  simulations.  The 
spatial  dimensions  of  the  root-cell  can  be  allowed  to  vary 
according  to  a  modified  constant-pressure  algorithm  due  to 
Berendsen,77  while  using  three-dimensional  pbcs  and  the  FMM 
algorithm.  Consequently,  bulk  samples  of  both  crystals  and 
dielectric  fluids  with  solvated  ionic  salts  at  any  pressure  or 
temperature  of  interest,  can  be  prepared  and  simulated. 

In  the  droplet  simulations,  the  basic  root-cell  of  the  FMM 
code  has  eight  tree  levels  with  the  smallest  leaf  cell  chosen  to 
be  1.2  nm  on  each  side.  This  gives  a  root-cell  of  dimension 
2(8-h  x  1.2  nm=  153.6  nm  dimensions  on  each  side.  In  studies 
pertaining  to  the  response  of  the  droplet  to  a  constant  external 
electric  field,  this  root-cell  was  replicated  three  times  along  the 
c-axis  (the  electric  held  direction)  to  give  a  total  tree-structure 
of  153.6  nm  x  153.6  nm  x  460.8  nm,  with  the  10  nm  diameter 
droplet  placed  in  the  center.  The  very  large  system  size  allows 
enough  space  for  the  parent  droplet  to  move  and  elongate  freely 
with  sufficient  time  for  charged  cluster  emission  to  occur.  Once 
the  simulation  begins,  no  temperature  control,  or  control  of  the 
droplet  center-of-mass  motion,  are  imposed, 

Given  the  size  of  the  present  droplets,  it  is  important  to  utilize 
an  efficient,  simple  representation  of  the  interatomic  forces, 
Formamide  is  treated  as  a  solid  body  with  quaternion  dynam¬ 
ics,77  implemented  by  using  a  midstep  implicit  leapfrog  algo¬ 
rithm78  that  has  been  shown  to  be  an  extremely  stable  integration 
scheme  with  minimal  energy  drifts  for  long  simulation  periods. 
The  geometry  of  the  formamide  molecule  is  taken  from  high- 
resolution  X-ray  studies  of  formamide  crystals.79  The  AMBER 
force  held  parameters80  were  used  to  describe  the  intermolecular 
van  der  Waals  interactions  between  formamide  molecules,  the 
interactions  between  the  formamide  molecules  and  the  atomic 
ions,  as  well  as  the  ion— ion  interactions.  CHELP-BOW  partial 
atomic  charges*1  were  used  for  the  overall  description  of  the 
electrostatic  potential  of  the  formamide  molecules.  Although 
long-range  Couloumbic  interactions  between  formamide  mol¬ 
ecules  were  not  truncated,  the  weak  van  der  Waals  interactions 
were  truncated  on  a  group  (molecular)  basis  through  the  use  of 
a  smooth  switching  function  described  by  Perera  et  al.82 

3.2.  Quantum  Chemical  Methods. 

3.2.1.  Solvated  Positive  Ions.  We  have  conducted  density 
functional  theory  (DFT)  calculations  of  the  observed  positive 
solvated-ion  structures  to  estimate  their  thermal  stability  and 
the  dissociation  energies  associated  with  loss  of  a  single 
formamide  solvent  molecule.  Calculations  were  performed  at 
the  B3LYP/6-3 1 0(d)  level  using  the  Gaussian  98/Gaussian  03 
quantum  chemistry  package.83  Geometries  were  optimized,  and 
normal-mode  frequencies  for  zero-point  energy  determinations 
were  calculated  for  solvated  ions  of  Na+(HCONH2)„  and 
K+(HCONH2)„,  where  n  =  1-6. 

3.2.2.  Solvated  Negative  Ions.  Initial  geometry  optimizations 
of  [l~(HCONH2)„],  n  —  1—5  were  performed  using  an  aug¬ 
mented  correlation  consistent  double-^  basis  set  for  second  row 
atom,84  6-3 1G  basis  set  for  hydrogen,85  Radom  ECP  basis  set 
for  iodine,86  and  the  B3LYP87  89  functional.  The  treatment  of 
negative  cluster  ions  required  a  higher  level  of  theory  due  to 
the  significantly  higher  degree  of  correlation,  and  thus  the 
geometries  were  refined  with  M0ller— Plesset  second-order 
perturbation  (MP2)  theory90  92  and  the  6-31  l+G(d,p)93  basis 
set  using  the  Radom  6-3 1  1  H-G(df)  basis  set  for  iodine.94 
Gradient  convergence  tolerance  for  geometry  optimizations  was 
set  to  lO'6  hartree/bohr.  Stationary  points  were  characterized 
by  a  positive  definite  Hessian  (matrix  of  energy  second 
derivatives).  Harmonic  zero-point  energy  and  temperature 
corrections  (to  298.15  K)  were  obtained  using  the  MP2  Hessians. 


All  calculations  were  performed  using  the  June  27,  2005  (R3) 
version  of  GAMESS  (General  Atomic  and  Molecular  Electronic 
Structure  System).9596 

3.3.  Unimolecular  Ion  Dissociation  Rate  Calculations.  The 
distributions  of  the  measured  ion  solvation  numbers  (that  is, 
the  distribution  function  of  the  number  of  solvating  forma¬ 
mide  molecules  in  a  solvated-ion  cluster)  deduced  from  the 
measurements,  are  compared  with  predictions  based  on  the 
assumption  that  the  ions  are  evaporated  from  a  liquid  at  or 
near  room  temperature.  The  internal  energy  of  a  solvated 
ion  increases  with  solvation  number,  and  the  dissociation 
energy  (that  is,  the  energy  associated  with  the  loss  of  a  single 
solvent  molecule)  declines  upon  increased  solvation.  Above 
a  certain  solvation  number,  the  average  internal  energy 
exceeds  the  dissociation  energy  of  the  solvated  ion,  leading 
to  high  dissociation  rates  in  relation  to  the  experimental  and 
simulation  time  scales.  Given  the  relatively  short  flight  times 
of  ^30  /.is  of  ions  in  our  mass  spectrometer,  it  is  possible 
that  metastable  ions,  i.e.,  those  with  internal  energies  above 
the  dissociation  limit,  reach  the  detector  prior  to  dissociating. 
To  assess  the  fractions  of  the  observed  ions  that  are  produced 
with  internal  energies  above  the  dissociation  limit,  and  to 
substantiate  the  present  interpretation  that  ion  held  evapora¬ 
tion  is  comparable  to  a  thermal  evaporative  process,  we  have 
conducted  RRKM  unimolecular  dissociation  rate  calculations 
of  ions  near  the  thermal  stability  limit  at  room  temperature. 

Using  the  normal-mode  frequencies  from  the  ab  initio 
calculations,  we  determine  the  thermal  internal  energy  distribu¬ 
tions  of  the  observed  solvated  ions  from  the  relation 

F(EiM)  d£mt  =  E,JEml)  Fvib(Evjb)  d£rot  dfvib  (8) 

where  Frol(£r()t)  and  FVib(£vib)  are  the  normalized  rotational  and 
vibrational  Maxwell— Boltzmann  energy  distribution  functions 
for  a  temperature,  T \  and 

£'im  =  Eroc  +  £vib  (9> 

For  those  ions  with  internal  energies  exceeding  the  binding 
energy,  we  determine  the  unimolecular  dissociation  rate  using 
a  RRKM  approach,  following  the  methodology  developed  by 
Hase  and  co-workers.97  Due  to  the  large  dipole  of  formamide 
(3.7  D)  and  the  resulting  strong  attractive  interaction  between 
the  solvated  ion  and  the  solvating  ligands,  there  is  no  centrifugal 
barrier  to  dissociation  and  a  very  loose  transition  state  is 
assumed.  We  assume  conditions  near  the  phase  space  limit,  and 
the  transitional  vibrational  modes  are  treated  as  rotors,  following 
Rodgers  et  al.98  All  systems  considered  include  a  single  solvation 
shell  and,  given  the  shallow  global  minima  and  near  symmetry 
of  the  normal-mode  frequencies,  all  ligands  are  treated  equiva¬ 
lently.  All  rotors  are  treated  as  active  in  both  the  reactant  and 
transition  states,  vibrational  states  are  treated  as  harmonic 
oscillators,  and  the  sum  and  density  of  states  are  determined 
using  the  Beyer—  Swinehart  direct  count  algorithm.99  Adiabatic 
treatment  of  both  reactant  and  transition-state  rotors  does  not 
markedly  change  the  results.  Furthermore,  given  the  “loose” 
nature  of  the  transition  state,  we  find  that  the  results  are  not 
sensitive  to  the  definition  of  the  transition  state  (i.e.,  the  distance 
of  the  dissociating  ligand  with  respect  to  the  ionic  center). 

The  survival  fraction,  P(n,T),  of  ions  of  a  particular  solvation 
number,  n,  at  a  temperature,  T,  is  then  determined  from  the 
expression 

P(n,T)  =  Ps(n,T)  +  f  PJn.T.Em) d Em  (10) 

where  Ps{n,T)  is  the  fraction  of  stable  ions  with  internal  energies 
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TABLE  1:  Electrospray  Source  Parameters  of  the 
Experimentally  Investigated  Nal— Formamide  Solutions0 


Nal  mole  ratio:  1:8.4  136.9 

conductivity:  1.64  S/m  1 .26  S/m 


extraction  voliage  (V) 

2000 

2000 

2400 

2400 

polarity 

+ 

- 

+ 

- 

/(total)  (nA) 

250 

464 

120 

200 

/(Faraday  cup)  (nA) 

28 

58 

63 

60 

m(QCM)(ng  s_l) 

2  4 

3.0 

5.4 

4.0 

w/r/(Daltons) 

8410 

4930 

8330 

6360 

£ma*(^//nm) 

1.70 

1.64 

1,42 

1.49 

^mai(nm) 

2.30 

2.48 

3.29 

2.98 

°  /(total)  and  /(Faraday  cup)  correspond  to  the  total  emission  and 
Faraday  cup  currents,  respectively.  m(QCM)  is  the  mass  flow  as 
measured  on  the  QCM.  E ma^  and  correspond  to  the  maximum 
electric  field  strength  and  the  associated  radius  of  curvature  at  the 
jet  neck  given  by  eqs  6  and  7. 

below  the  dissociation  limit,  Pm{nyT,Eml)  d£mt  is  the  fraction 
of  ions  with  energies  in  an  interval  Eint  -F  d£m{;  k{Emt)  is  the 
ummolecular  dissociation  rate  constant,  and  r  is  the  average 
flight  time  of  the  respective  ions  in  the  mass  spectrometric 
apparatus  described  in  section  2.  The  experimental  flight  times 
are  on  the  order  of  25  //s. 

4.  Results 

4.1.  Mass  Spectrometric  Measurements.  The  concentrated 
and  dilute  Nal  — formamide  solutions  were  electrosprayed  tn  a 
Taylor  cone-jet  mode,  The  operational  parameters  of  the 
experiments  are  summarized  in  Table  1.  For  the  1:8.4  mol  ratio 
solution,  the  extraction  voltage  was  2000  V  for  both  polarities, 
whereas  a  slightly  higher  voltage  was  necessary  for  stable 
operation  of  the  more  dilute  1:36.9  solution.  The  measured 
charge  emission  currents  correspond  to  /(total)  in  the  Table  1. 
Mass  flows  registered  by  the  QCM  detector  are  also  listed.  The 
average  charged  particle  mass-to-charge  ratio,  mlq ,  is  then 
derived  from  the  ratio  of  mass  flow  to  Faraday  cup  current. 
The  comparison  between  the  total  emission  current  and  the 
Faraday  cup  current  which  collects  at  a  solid  angle  close  to 
that  of  the  QCM,  allows  an  assessment  of  the  electrospray 
divergence.  The  markedly  higher  divergence  observed  for  the 
concentrated  solution  can  be  attributed  to  the  lower  acceleration 
voltage,  the  higher  current  density,  as  well  as  the  possibility  of 
changing  electrospray  dynamics.  This  will  be  further  addressed 
in  the  discussion  (see  section  5). 

The  relatively  low  mlq  values  suggest  that  both  solutions  are 
sprayed  in  a  mixed  droplet/ion  emission  regime,  where  the  main 
fraction  of  mass  flow  is  carried  by  the  droplets,  and  a  substantial 
fraction  of  the  current  comes  from  tons.  This  indicates  that  the 
present  Taylor  cone-jets  operate  at  dimensions  where  ion  field 
evaporation  occurs.  Field  evaporation  is  further  verified  by  the 
determination  of  the  maximum  electric  field  strength  and  the 
corresponding  characteristic  radius  of  curvature  according  to 
eqs  6  and  7.  The  flow  rates,  Q ,  were  obtained  from  the  mass 
flow  assuming  a  density  of  1.13  g  cm-3.  The  surface  tension  of 
pure  formamide,  0.0575  N  m  *,  was  assumed  for  the  calculation 
of  Note  that  eqs  6  and  7  are  based  on  the  total  flow  rate 
so  that  the  derived  values  based  on  the  QCM  registered  mass 
flow  must  be  regarded  as  upper  (£ma*)  and  lower  ( Rmax )  limits. 
Depending  on  the  charge  emission  dynamics,  it  is  possible  that 
nonuniformity  exists  in  the  angular  mlq  distributions,  thereby 
causing  differences  hetween  mlq  values  determined  on  axis,  as 
with  the  present  near-field  probe,  compared  to  those  ohtained 
from  the  total  current  and  mass  flow.  The  representative 
maximum  electric  field  strength  ranges  from  1.42  to  1.7  V/nm, 


Mass  (amu) 

Figure  4.  Sample  positive  and  negalive  ion  mass  spectra  observed 
from  the  vacuum  electrospray  of  a  concentrated  (28.4  wt  %) 
Nal— formamide  solution  Solvation  number  assignment  are  shown 
for  solvated  and  complex  ion  structures. 

with  higher  fields  derived  for  the  concentrated  solution,  whereas 
the  jet  dimensions  are  associated  with  radii  in  the  range  between 
2.3  and  3.29  nm,  with  the  largest  radii  determined  for  the  dilute 
solutions. 

The  negative  total  emission  currents  tend  to  be  almost  twice 
as  high  as  those  observed  for  positive  charge  emission  Two 
significant  differences  are  observed  for  the  dilute  solution  in 
comparison  with  the  concentrated  solution:  (i)  the  total  emission 
current  is  lower,  and  the  fractional  current  registered  on  the 
Faraday  cup  is  higher,  signifying  a  less  divergent  spray,  and 
(ii)  the  negative  ion  average  mass-to-charge  ratio  is  higher  for 
the  dilute  solution,  whereas  the  positive  ion  mass-to-charge  ratio 
does  not  change  significantly  upon  dilution. 

Figure  4  shows  positive  and  negative  ion  mass  spectra 
recorded  for  the  1:8.4  Nal— formamide  solution  The  spectra 
have  been  corrected  w  ith  respect  to  the  mass  spectral  transmis¬ 
sion  function.  The  low  mass  resolution  can  be  attributed  to  the 
high  transmission  energies  of  550  eV.  The  most  abundant 
positive  ion  mass  can  be  attributed  primarily  to  solvated  sodium 
ions,  Na+(HCONH2)„,  with  the  n  =  4  solvated  ion  being  the 
most  abundant  positive  ion  species.  Trace  abundances  of  ions 
attributed  to  solvated  ions  including  an  ion  pair, 
Na+(Nal)(HCONH2)„,  are  also  observed.  Such  complex  ions 
are  not  detected  with  a  positive  polarity  in  the  electrospray  mass 
spectrum  of  the  1:36.9  solution.  The  mass  resolution  is 
insufficient  to  fully  resolve  the  ions  including  an  ion-pair  from 
the  simple  solvated  ions.  The  individual  abundances  are 
determined  by  deconvoluting  the  data  using  fits  of  Gaussian 
peak  shapes  given  by  the  experimental  resolution  and  centered 
at  the  expected  mass  positions.  The  derived  normalized  abun¬ 
dances  are  shown  for  the  two  concentrations  in  Figure  5a,b. 
The  abundances  of  the  simple  solvated  ions  are  essentially 
identical  for  the  two  concentrations  except  for  higher  abundances 
at  the  highest  solvation  numbers  in  case  of  the  more  dilute 
solution. 

When  operated  in  a  negative  mode,  the  main  ohserved  ions 
are  solvated  iodine  anions,  1  (HCONH2)*,  with  a  distribution 
peaking  at  n  =  2  in  the  mass  spectrum  observed  for  the 
concentrated  (1:8.4)  solution  (Figure  4).  Solvated  ions  including 
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Solvation  Number,  n 

Figure  5,  Na+  solvation  number  distributions  derived  from  the  mass  spectra  of  (a)  1:8.4  (top)  and  (b)  1 :36.9  (bottom)  mole-ratio  Nal- forma  nude 
solutions.  I  solvation  number  distributions  derived  from  the  mass  spectra  of  (c)  1:8.4  (top)  and  (d)  1  36.9  (bottom)  mole-ratio  Nal  forma m id e 
solutions. 


ion-pairs  1  (Nal)(HCONH2)„,  are  also  observed.  The  abundances 
following  deconvolution  are  shown  in  Figure  5c, d.  Contrary  to 
the  Na+  positive  tons,  dilution  has  a  marked  effect  on  the 
observed  solvated  1“  mass  spectrum.  The  distribution  of  simple 
solvated  tons  becomes  significantly  broader  upon  dilution,  and 
it  peaks  at  a  higher  solvation  number,  n  =  5.  As  with  the  positive 
ions,  the  complex  ions,  l”(Nal)(HCONH2)„,  are  not  detected 
upon  dilution. 

Although  no  simulations  were  carried  out  for  K1  solutions, 
we  found  it  useful  for  comparative  purposes  to  also  investigate 
mass  spectrometrically  Kl— formamide  solutions.  The  Kl  — 
formamide  solution  could  not  be  electrosprayed  at  the  same 
concentration  (28.4  wt  %)  used  for  the  Nal— formamide  solution 
due  to  problems  caused  by  crystallization  in  the  capillary  orifice 
region.  Stable  Tay  lor  cone -jets  were  produced  with  Kl-forma- 
mide  mole-ratios  of  1  ;41  and  1 : 83  at  extraction  voltages  of  2100 
and  2300  V,  respectively.  Simple  and  complex  solvated  tons 
were  observed,  as  for  the  Nal -formamide  system.  The  a  bun 
dances  are  shown  in  Figure  6a— d.  Similar  to  the  Na+  solvated 
ion  distributions,  dilution  results  in  only  a  minor  change  in  the 
K+  solvation  numher  distributions.  A  small  shift  toward  lower 
solvation  numbers  is  found  upon  dilution.  However,  the 
distributions  (that  peak  near  n  =  5)  extend  to  significantly  higher 
solvation  numbers  compared  to  the  Na+  distributions  discussed 
above,  and  they  may  continue  to  even  larger  solvation  numhers 
but  are  truncated  by  the  range  of  our  quadrupole  mass 
spectrometer.  Complex  ions,  K+(KI)(HCONH2),„  are  generally 
more  ahundant  than  in  the  Nal  solutions,  and  they  appear  to  he 
more  abundant  for  the  dilute  solution.  The  solvated  anion 
distributions  show  the  same  trend  of  peaking  at  higher  solvation 
numbers  upon  dilution  as  observed  for  the  Nal  solutions. 
Contrary  to  the  dilute  Nal  solution,  for  the  Kl  case  complex 
anions  appear  to  be  more  abundant  upon  dilution. 

Figure  7  compares  the  retarding  potential  measurements  of 
energy  distributions  per  unit  charge  of  mass-selected 
Na+(HCONH2)4  ions  with  high  mass  charges  ( mlq  >  400  amu). 


obtained  by  operating  the  quadrupole  in  a  radio-frequeney-only 
high-pass  mode.  I  he  capillary  potential  was  hiased  at  550  V 
while  the  retarding  grid  analyzer  was  at  ground.  ITie  data  show 
that  high  mass  charged  particles  correspond  most  likely  to 
charged  droplets  and  are  formed  at  a  potential  ^200  V  below 
the  capillary  or  bulk  liquid  potential,  whereas  the  solvated  ions 
are  formed  with  a  broader  energy  distribution  centered  slightly 
below  the  droplet  appearance  energy,  w  ith  a  minor  contribution 
at  higher  energies.  Very  similar  distributions  were  obtained  for 
ions  with  other  solvation  numbers;  at  the  same  time,  the  energy 
distributions  do  not  change  noticeably  for  negative  ions  and 
for  varying  dilutions. 

4.2.  Solvated  Ion  Thermochemistry  and  Ion  Survival 
Fractions.  Calculated  global  minimum  structures  of  several 
solvated  ions  are  shown  in  Figures  8  and  9;  the  geometrical 
parameters  for  the  optimized  structures  of  Na(HCONH2),/ ,  n 
=  1  —4,  and  l(HCONH2)„  ,  n  —  1  —5,  are  given  in  the  Appendix, 
Tahles  Al  and  A2,  respectively.  The  binding  energies  of  the 
calculated  solvated  ions,  Na(HCONH2)/,  n  =  1—6, 

K(HCONH2)„\  n  =  1-6,  and  l(HCONH2)„  ,  n  =  1-5, 
determined  from  the  quantum  chemical  calculations,  are  shown 
in  Table  2.  Also  listed  are  the  zero-point  energies.  E7V,  and  the 
internal  energies,  Emh  above  the  zero-point  energies  at  298. 1 5 
K  (difference  between  thermal  energy  and  zero-point  energy). 
It  is  seen  that  the  internal  energy  at  298.15  K  exceeds  the 
binding  energy  for  ions  with  a  solvation  number  above  n  =  4 
for  the  cationic  solvated  ions  (Na+  and  K-1-),  and  above  n  =  3 
for  the  1“  solvated  ions,  respectively.  Thus  solvated  ions  with 
n  above  the  solvation  numbers  limits  are  thermally  unstable  and 
given  sufficient  time  they  will  dissociate,  desorbing  a  formamide 
molecule. 

Lifetimes  of  ions  with  solvation  numbers  near  the  thermal 
stability  limit  were  calculated  for  Na+  and  1  solvated  ions. 
Table  3  lists  the  thermal  ion  survival  fractions  calculated  for 
selected  Na+(HCONH2)»  ions  evaporated  with  characteristic 
temperatures  of  300,  350  and  400  K.  The  survival  fractions  are 
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Figure  6.  K*  solvation  number  distributions  derived  from  the  mass  spectra  of  (a)  1:41  (top)  and  (h)  1:83  (hottoni)  mole-ratio  Kl  —  formanude 
solutions.  !  solvation  number  distributions  derived  from  the  mass  spectra  of  (c)  1 :41  (top)  and  (d)  1 :83  (bottom)  mole-ratio  Kl  —  forma niide  solutions. 
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Figure  7.  Retarding  potential  measurements  of  solvated  ion  and  droplet 
(corresponding  to  high  pass,  mlq  >  400)  energy  distributions.  The 
capillary  potential  was  biased  at  550  V,  and  the  analyzer  was  at  ground 

calculated  for  infinite  observation  times  and  mass  spectrometer 
flight  times  and  are  based  on  the  unimolecular  dissociation  rate 
calculations  described  in  section  3.3.  Tabic  4  lists  corresponding 
survival  fractions  for  selected  I  (HCONFb),,  ions. 

4.3.  Molecular  Dynamics  Simulations.  Neutral  dielectric 
and  charged  droplets  and  their  response  to  uniform  electric  fields 
have  been  of  great  interest  for  a  rather  long  time  for  both 
fundamental  scientific  reasons  and  technological  ones. 8,17  71-74 
The  salient  features  of  the  physics  associated  with  electrosprays 
operating  in  the  cone-jet  mode  may  be  seen  already  in  the 
simpler  problem  of  a  conducting  liquid  droplet  elongated  tn  a 
uniform  external  electric  field,  including  the  formation  of 
electrified  jets,  as  well  as  the  emission  and  acceleration  of 
charged  species  at  the  ends  of  the  elongated  parent  drops. 
Indeed,  Taylor17  studied  the  instabilities  and  jet  formation  of 
fluids  emerging  from  an  electrospray  device  as  a  way  to 
investigate  in  a  more  controlled  manner  the  physics  related  to 
the  elongation  and  jet  formation  of  droplets  in  a  uniform  electric 
field. 


In  light  of  the  above  we  initiated  atomistic  molecular 
dynamics  simulation  studies,  aiming  at  exploration  of  held- 
induced  fluid  instabilities  relevant  to  the  operation  of  electro¬ 
spray  devices  hy  studying  first  the  closely  related  phenomena 
present  in  the  behavior  of  dielectric  droplets,  with  and  without 
dissolved  salts.  As  discussed  below,  the  results  of  the  MD 
simulations,  which  are  found  to  agree  well  with  the  experimental 
electrospray  results  measurements,  provide  deep  insights  into 
the  physical  origins  of  the  ohserved  behavior  As  a  prelude  to 
our  discussion  of  the  Nal— formamide  solution  droplets  that  are 
our  primary  focus,  we  present  in  section  4.3.1  results  for  pure 
formamide  droplets  which  are  taken  to  have  the  same  size  as 
the  Nal— formamide  solution  droplets  which  we  analyze  in 
section  4.3.2. 

4.3.1.  Simulations  of  Pure  Formamide.  Because  the  inter¬ 
atomic  and  intennolecular  interaction  potential  parameters 
employed  in  our  simulations  were  not  specifically  developed 
to  reproduce  certain  physical  properties  of  the  simulated 
formamide  (such  as  the  melting  point  and  dielectric  constant), 
we  performed  first  exploratory  simulations  that  allow  us  to  assess 
the  physical  properties  of  the  simulated  liquid  and  the  faithful¬ 
ness  of  the  simulation  model.  Constant  pressure  simulations  and 
liquid— solid  coexistence  studies100  indicate  that  our  description 
of  formamide  yields  a  liquid  density  that  is  to  within  a  few 
percent  of  the  experimental  value,  and  a  melting  point  that  is 
approximately  10  K  above  the  experimental  value  of  275  K. 
For  an  improved  comparison  with  experiment  and  to  avoid 
freezing  of  the  field-evaporating  salt  solutton  droplets,  the 
simulations  are,  therefore,  performed  at  310  K  (that  is,  10  K 
above  the  experimental  conditions). 

The  droplets  used  in  our  simulations  were  created  hy  “carving 
out”  a  spherical  region  from  an  equilibrated  bulk  liquid 
simulation  that  was  performed  at  constant  temperature  and  zero 
pressure,  with  a  subsequent  equilihration  prior  to  the  application 
of  an  external  field  The  diameter  of  the  spherical  droplet  chosen 
for  the  MD  simulations  is  10  nm  and  it  contains  7150  formamide 
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Figure  8.  (a)-(d)  Na^HCONHi)*,  n  =  1-4,  optimized  solvated  ion 
structures  calculated  at  the  B3LYP/6-31G(d)  level. 


molecules  (each  molecule  composed  of  six  atomic  partial 
charges).  This  droplet  is  large  enough  to  allow  meaningful 
comparison  with  continuum  electrohydrodynamic  theories,  and 
in  addition  its  curvature  is  similar  to  the  estimated  local 
curvature  of  the  experimentally  probed  liquid  jet  at  the  location 
where  the  electric  held  strength  is  maximal  (see  eqs  6  and  7). 
From  a  practical  point  of  view,  the  chosen  droplet  size  is 
sufficiently  small  to  allow  a  wide  range  of  computer  experiments 
to  be  performed  in  a  reasonable  amount  of  time. 

The  droplet  was  equilibrated  in  a  range  of  increasing  electric 
field  strengths  (from  0  to  3  V/nm).  The  effect  of  the  externally 
imposed  electric  field  is  depicted  in  Figure  10,  where  we  display 
snapshots  of  the  droplet  for  three  values  of  the  applied  field, 
illustrating  instability  at  a  critical  electric  field  strength.  At  field 
strengths  below  ^0.5  V/nm,  the  droplet  undergoes  mild 
elongation  with  increasing  field  strength.  At  a  field  strength  of 
0.625  V/nm  the  droplet  becomes  unstable  and  undergoes  a 
pronounced  elongation,  resulting  in  a  change  in  aspect  ratio  // 


Figure  9.  (a)— (e)  1  (HCONH>)„,  /i  =  1—5,  optimized  solvated  ion 
structures  calculated  at  ihe  MP2/6-31  l+G(d,p)  level. 


=  c/a  (where  c  and  a  are  the  semimajor  and  semi  minor  axes  of 
a  prolate  ellipsoidal  fit  to  the  droplet  shape)  from  //  ~~  2  to  //  ~ 
12.  The  elongation  is  also  reflected  in  changes  in  the  average 
molecular  dipole  moment  projected  onto  the  z  axis.  There  have 
been  a  number  of  studies  addressing  the  critical  external  field 
at  which  a  dielectric  droplet  displays  this  instability.  In  early 
work,  Taylor17  theoretically  predicted  that  a  droplet  of  initial 
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TABLE  2:  Thermochemistry  of  Na^(HCONH2)„,  n  =  1—6, 


K+(HCONH2)„,  n  =  1  6,  and  I  (HCONH2)„,  n  =  1-5, 
Solvated  Ions  Determined  with  the  Indicated  ab  Initio 
Methods" 

B3LYP/6-31G(d) 

n 

E/p c 

Em 

D(MF„  ,-F) 

Na+(HCONHi),. 

1 

30.1 

3.33 

39.7 

2 

59.9 

6.71 

33.3 

3 

89.4 

10.3 

24.1 

4 

1 18.4 

14.0 

17.2 

5 

149.5 

16.3 

15.7 

6 

178.2 

20.2 

7.2 

B3LYP/6-3lG(d) 

n 

E,pc 

E\n\ 

D(MF„  ,-F) 

K+(HCONH2)„ 

1 

29  8 

3.46 

28.7 

59.3 

7.02 

24.0 

3 

88.7 

10.6 

19.5 

4 

117.7 

14.4 

15.6 

5 

148.8 

16.8 

14.4 

6 

178.3 

20.3 

10.8 

MP2/6-3 1 1  +G(d,p) 

n 

l'/pc 

Fint 

D(MF„  ,-F) 

1 

27.8 

1  (HCONH2)„ 

3.6 

17.4 

2 

56.0 

7.0 

14.9 

3 

83.7 

10.8 

12  1 

4 

111.8 

14.3 

10.9 

5 

140.0 

17.1 

8.4 

"  The  internal  energy.  Eml.  signifies  the  298.15  K  thermal  energy 
above  the  /ero-point  energy,  E/pe.  The  binding  energy, 
IX MF„  i  —  F),  refers  to  the  loss  of  a  single  solveni  molecule. 
F^HCONFL.  All  values  are  given  in  kcal/mol. 


TABLE  3:  Survival  Fractions  for  Selected  Na^(HCONH2)„ 
Ions  Evaporated  with  Characteristic  Temperatures  of  300, 
350  and  400  K" 


t 

n  =  4 

n  =  5 

ti  —  6 

T  =  300  K 

CO 

0.927 

0.503 

8.0  x  10  5 

MS  (26,  28,  31  p s) 

0.994 

0.662 

1  0  x  10  4 

T  =  350  K 

oo 

0.702 

0.148 

0.000 

MS  (26,  28,  31  [is ) 

0.920 

0.279 

0.000 

T  =  400  K 

oo 

0.332 

0.035 

0.000 

MS  (26,  28,  31  ps) 

0.718 

0.070 

0.000 

“  The  survival  fractions  are  listed  for  infinite  observation  times,  f, 
and  flight  limes  through  mass  specirometrie  experiment  (MS). 


radius  R  (in  zero  external  field)  becomes  unstable  at  an  external 
field  strength,  Ecrn\c* i,  satisfying  (MKS  units): 

^critical  =  '  .625(y/2;7T0c/)12  (11) 

Using  y  =  0.046  N/m  for  the  surface  tension  of  the  simulated 
droplet  (computed  in  separate  simulations)100  and  d  =  9.68  nm 
for  the  computed  equimolar  diameter  results  in  a  predicted 
cntical  applied  electric  field  E  =  0.48  V/nm,  which  is  only 
slightly  lower  than  the  interval  of  0.5-0.625  V/nm  in  which 
the  instability  is  observed  to  occur  in  the  MD  simulation;  we 
note  that  if  the  experimental  surface  tension101  y  =  0.0575  N/m 


TABLE  4:  Survival  Fractions  for  Selected  I  (HCONH2)„ 
Ions  Evaporated  with  Characteristic  Temperatures  of  300, 
350  and  400  K" 


r 

n  =  3 

n  =  4 

=  5 

oo 

MS  (29,  32,  34  /( s) 

r  =  300  K 
0.859 

0.957 

0.295 

0548 

0.036 
0.09 1 

oo 

MS  (29.  32,  34  /<s) 

T  =  350  K 
0.620 

0.811 

0.088 

0.224 

0.006 

0.016 

OO 

MS  (29,  32,  34  ps) 

T  =  400  K 
0.341 

0.570 

0.023 

0.067 

0.000 

0.001 

"  The  survival  fractions  are  listed  for  infinite  observation  times,  t , 
and  flight  times  through  mass  spectrometne  experiment  (MS). 
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Figure  10.  Atomic  configurations  laken  from  a  molecular  dynamics 
simulaiion  of  a  10  nm  formamide  droplei  prior  lo  and  afier  the 
application  of  an  electric  field.  The  formamide  molecules  are  depicted 
in  blue  When  a  crilical  field  is  exceeded,  the  droplet  deforms  lo  an 
ellipsoidal  stable  configuration  The  crilical  field  is  calculated  for  a  10 
nm  diameier  droplei  and  assumes  a  surface  tension  of  0.046  N/m.  For 
larger  fields,  highly  elongated  shapes  with  sharply  pointed  ends  develop 
along  ihe  field  direction  The  aspect  raiios,  >/,  for  the  various 
configurations  are  specified.  The  Z-axis  points  in  the  electric  field 
direction. 


is  used,  the  predicted  critical  field  strength  is  0.53  V/nm.  We 
conclude  that  a  continuum  formalism  provides  an  adequate 
description  of  our  results  obtained  for  nanometer  scale  droplets. 
Additionally,  it  was  discovered010  that  at  higher  electric  fields 
(about  1.4  V/nm)  field-induced  crystallization  of  the  pure 
formamide  droplet  occurs. 

4.3.2.  Simulations  of  10  nm  Naf  —formamide  droplets. 
Having  illustrated  that  the  thermodynamics  and  phase  behavior 
of  a  dielectric  droplet  can  be  modified  and  manipulated  through 
the  use  of  external  electric  fields,  we  turn  next  to  investigations 
of  formamide  droplets  with  a  dissolved  Nal  salt  at  concentrations 
that  are  comparable  to  those  used  in  the  experiments  described 
in  section  4  1.  Simulations  were  performed  on  10  nm 
Nal— formamide  solutions  with  Nal— formamide  ratios  of  1:8 
and  1:16,  respectively  in  fields  of  0.5,  0.625,  0.75,  1.0  and  1.5 
V/nm  The  concentrated  droplets  contain  7150  formamide 
molecules  and  900  dissolved  Nal  units,  corresponding  to  a  29 
wt  %  solution  that  is  nearly  identical  to  the  concentrated  solution 
used  in  the  experiments.  The  salt  solution  droplets  were  prepared 
similarly  to  the  pure  formamide  droplets  described  above. 

Snapshots  of  the  1:16  salt  solution  droplet  subjected  to  an 
electric  field  of  0.625  V/nm  are  shown  in  Figure  1 1  A  stage  of 
initial  elongation  along  the  field  direction  is  shown  in  Figure 
1  lb,  400  ps  after  the  field  has  been  turned  on.  Figure  1  lc,  taken 
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before  field  is  applied  (0  ps)  initial  elongation  (400  ps) 


Figure  11.  Atomic  configurations  taken  from  a  molecular  dynamics 
simulation  of  a  10  nm  diameter  droplet  of  a  I  16  Nal  — formamide 
solution,  selected  during  different  stages  of  elongation  following  the 
application  of  a  uniform  external  electric  field  of  0.625  V/nm.  Na+ 
and  I  ions  are  shown  as  red  and  yellow'  spheres  respectively,  and 
formamide  (F)  molecules  are  in  blue.  In  the  inset  for  //  =  2,  showing 
a  hydrogen-bond  connected  “molecular  wire”,  the  oxygen  atom  of  the 
formamide  molecule  is  shown  in  pink,  the  hydrogens  in  light  blue,  the 
carbon  in  light  green,  the  nitrogen  in  dark  blue,  (d)  and  (e)  show 
formation  of  a  nanojet  and  subsequent  emission  from  its  tip  of  nanoscale 
clusters  of  ions  solvated  by  formamide  molecules  Time  is  measured 
(in  ps)  with  respect  to  the  instant  when  the  field  is  turned  on  (/  =  0): 
(a)  0:  (b)  4(X)  ps;  (e)  700  ps:  (d)  1 100  ps,  (e)  1125  ps  The  aspect 
ratios,  //,  for  the  various  configurations  are  also  indicated. 


300  ps  later,  exhibits  further  elongation  of  the  droplet,  and  the 
onset  of  charge  emission  in  the  form  of  a  Na+  solvated  ion.  In 
the  snapshot  shown  in  Figure  1  Id,  which  was  recorded  1.1  ns 
after  the  initial  application  of  the  field,  formation  of  jets 
emanating  from  both  ends  of  the  elongated  droplet  are  apparent. 
Although  a  large  solvated  ion  (with  ~20  formamide  ligands) 
with  a  single  Na+  ion  at  the  center,  is  seen  to  be  on  the  verge 
of  being  emitted  at  the  positive  end  of  the  parent  droplet  (on 
the  right  in  Figure  lid),  a  large  complex  negative  ton, 
1  (NalKHCONFb)/;,  is  seen  to  be  detaching  at  the  other 
(negative)  end  of  the  droplet  (left  in  Figure  1  Id).  About  25  ps 
after  the  instant  when  the  configuration  shown  in  Figure  1  Id 
was  recorded,  a  section  of  the  material  streaming  from  the 
positively  charged  end  of  the  droplet  thins  down  to  a  narrow 
hydrogen  bonded  filament  of  formamide  molecules.  This  results 
in  the  detachment  of  a  cluster,  or  nanodroplet,  of  charge  4-3 
(see  Figure  lie)  consisting  of  17  Na+,  14  1“  and  ~200 
formamide  molecules.  This  corresponds  to  a  mass-to-charge 
ratio  of  about  1  1  2  kamu  (kDaltons). 

The  charged  clusters  that  detach  from  both  ends  of  the  droplet 
are  accelerated  by  the  applied  electric  field.  Their  physical 
characteristics  (charge,  mass,  size,  energies)  are  recorded  when 
they  cross  two  planes  (which  we  refer  to  as  “detection  planes”) 
located  at  ±200  nm  from  the  droplet  center  along  the  field 
direction  (the  z-axis).  Properties  of  the  recorded  positively  and 
negatively  charged  clusters  are  summarized  in  Figure  12  as  a 
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Figure  12.  Properties  of  the  charged  clusters  emitted  from  the  droplet 
described  in  Figure  I  I  as  a  function  of  cluster  radius:  (a)  the  estimated 
electric  field  at  the  droplet  surface,  (b)  the  charge-to-mass  ratios  and 
(c)  the  observed  (MD)  charge  divided  by  the  predicted  charge  for 
Rayleigh  instability  and  fission.  The  dashed  and  solid  lines  compare 
MD  results  to  estimates  based  on  the  experimental  clectrospray  results 
(from  Table  I;  see  text  for  details). 


function  of  the  cluster  radii,  /?,  estimated  by  using  the  equimolar 
radius  defined  by  the  cluster  mass,  m,  the  computed  density  of 
the  Nal— formamide  solution,  p  ~  1.32  g/cm\  and  the  relation 
m  —  4/3.T 7?3p.  Figure  12a  displays  the  estimated  surface  electric 
field,  E,  of  each  cluster  using  Coulomb's  law  (eq  2).  A  change 
in  the  dependence  of  E  on  the  cluster  radius,  R,  is  observed  at 
R  %  I  nm,  with  clusters  having  radii  R  <  1  nm  corresponding 
to  single  ions  solvated  by  formamide  molecules,  whereas  clusters 
with  R  >  1  nm  are  characterized  as  multi-ion  droplets.  In  the 
latter,  dispersion  forces  are  expected  to  have  a  more  important 
contribution  in  conjunction  with  electrostatic  interactions.  The 
change  in  behavior  occurring  for  clusters  with  R  >  1  nm  is 
reflected  also  in  the  qlm  and  17/^Rayicigh  ratios  plotted  in  Figure 
12b,c,  respectively.  The  arrow  on  the  horizontal  axis  in  Figure 
12  corresponds  to  the  experimental  estimate  (see  Table  1  for 
the  1 :8.4  mol  ratio)  for  the  radius  of  curvature  of  the  electrospray 
jet,  /?mux,  at  the  point  of  maximal  electric  field,  EmdX  (see  eqs  6 
and  7).  Essentially,  all  the  clusters  recorded  in  the  MD 
simulation  have  radii  smaller  than  Emax.  The  few  clusters  with 
radii  larger  than  7?milx  appeared  at  the  very  end  of  the  simulation 
when  breakup  of  the  droplet  into  large  fragments  was  observed. 
The  experimental  estimate  for  Em ax  =  1.67  V/nm  (determined 
from  eqs  6  and  7  as  the  average  between  the  positive  and 
negative  polarities)  is  indicated  by  the  solid  line  in  Figure  12a. 
The  estimated  surface  field  of  the  simulated  clusters  remains 
very  close  to  Emax  for  cluster  radii  in  the  droplet  regime. 

Figure  12b  shows  the  charge-to-mass  ratios  for  the  clusters 
emitted  by  the  simulated  droplet.  The  mean  value  of  the 
experimental  qlm  (Table  1 ),  averaged  over  polarities,  is  shown 
as  the  horizontal  dashed  line  (near  the  bottom  in  Figure  12b). 
The  MD  simulation  values  rapidly  approach  the  mean  experi¬ 
mental  value  as  the  cluster  sizes  enter  the  droplet  regime  If 
one  uses  £max  as  an  estimate  for  the  limiting  value  of  the  surface 
field  for  a  charged  droplet  in  eq  2  (Coulomb  law)  to  write  q  = 
4,T/?2foEmax,  an  estimate  for  the  upper  bound  of  q/m  can  be 
obtained  from 
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q/m  =  3F0EmJpR  (12) 

The  upper  bound  for  qlm  as  a  function  of  Emax  and  R  is  indicated 
in  Figure  12b  by  a  solid  line.  It  is  evident  that  the  data  for 
charged  droplets  (R  >  I  nm)  recorded  in  the  MD  simulations 
follow  the  predicted  trend 

In  Figure  1 2c  we  display  the  ratio  between  the  observed  (MD) 
charge  and  the  predicted  charge  for  Rayleigh  instability  and 
fission,  computed  from  eq  1 ,  with  the  surface  tension  y  =  0.046 
N/m  determined  from  the  MD  pure  formamide  droplet.  The 
charge  for  clusters  in  the  droplet  regime  exhibits  variations  about 
0  5  of  the  Rayleigh  limit.  Using  eqs  1  and  2,  one  obtains  the 
upper  limit  for  the  ratio,  ^Rayleigh- 

^Rayle.gh  =  ('„ (13) 

The  corresponding  values  of  eq  13  are  shown  as  a  solid  line  in 
Figure  12c.  This  simple  estimate  captures  the  general  trend  and 
values  observed  in  the  MD  simulation.  As  mentioned  earlier, 
the  largest  droplets  were  produced  at  the  end  of  the  simulation 
during  the  disintegration  of  the  main  droplet. 

Profiles  of  ion  (Na+,  solid  line;  I  ,  dashed  line)  and  solvent 
(circles)  molecular  densities  along  the  semimajor  axis  of  the 
droplet  spheroids  for  the  cases  (with  rj  =  1,  1.5  and  2)  shown 
in  Figures  1  la— c,  are  compared  in  Figure  13.  Also  included  in 
this  figure  are  profiles  of  the  local  electric  field  (squares  for  the 
field  outside  the  droplet  and  dotted  line  for  the  interior  field) 
and  the  ionic  charge  density,  the  latter  given  by  the  difference 
between  positive  and  negative  ion  densities  (see  the  lower  solid 
line,  near  the  horizontal  axis,  in  each  panel).  The  plotted  electric 
field  does  not  include  the  external  field  (0.625  V/nm).  The 
densities  are  calculated  by  defining  prolate  spheroid  shells  that 
match  the  instantaneous  droplet  shape  and  end  curvatures,  and 
by  counting  the  respective  species  between  neighboring  shells 
and  within  a  1  nm  core  radius  with  respect  to  the  c-axis. 
Although  the  radial  diameters  of  the  narrowing  droplet  ends 
(along  the  e-axis)  are  less  than  that  of  the  sampled  core  diameter 
(2  nm),  the  contributions  to  the  statistics  of  the  density  come 
only  from  material  inside  the  prolate  spheroid  shells,  thereby 
yielding  an  accuracy  superior  to  that  obtained  via  simple  binning 
of  the  c-axis.  This  process  was  not  implemented  for  the  case  of 
the  droplets  in  Figures  1  ld,c  due  to  the  narrow  ,  rapidly  vary  ing, 
curvature.  The  data  are  obtained  from  300  ps  time  averages  in 
the  vicinity  of  the  snapshots  shown  in  Figures  1  la— c. 

The  electric  field  strength  along  the  z-ax is  is  computed  by 
binning  (during  the  simulation)  the  number  of  each  species  (Na, 
I,  and  the  six  atoms  in  each  formamide  molecule,  each  having 
their  own  partial  charge)  with  the  z-  and  r-axis  having  bin  widths 
of  0.2  nm.  The  binned  data  is  averaged  over  a  specified  time  in 
which  the  droplet  does  not  change  its  shape  or  size  significantly. 
On  the  basis  of  the  charge  of  each  species,  the  binned  number 
of  atoms  provide  the  average  amount  of  charge  coming  from 
each  species  in  each  of  the  radial  “rings"  associated  with  each 
z-bin.  The  potential  V(Z)  at  a  point,  Z,  on  the  e-axis  is  given  as 
a  sum  of  the  contributions  from  the  total  charge  qu  in  each  of 
the  binning  rings  of  radius,  r\,  for  all  the  e-bins  (except  for  the 
central  radial  bin  in  the  e-bin  for  which  the  field  is  determined): 


The  electric  field  at  point  Z,  originating  from  the  ionic  and 
molecular  charges  (not  including  the  external  field)  is  then  given 
by  E{Z)  =  -d  V{Z)!d  Z, 

The  field-free  droplet  (Figure  13a)  exhibits  a  depletion  of 
ions  at  the  interface,  with  a  slightly  higher  ion  density  in  the 


Figure  13,  Number  densities  and  induced  eleetrie  field  of  lhe  1:16 
ratio  Nal— formamide  droplet  along  the  field  direction  for  (a)  no  external 
field,  (b)  elongaiing  in  an  applied  (external)  field  of  0.625  V/nm,  before 
emitting  ions,  and  (e)  during  the  time  of  initial  cluster-ion  emissions. 
The  profiles  correspond  to  lhe  respective  snapshots  in  Figure  1  la— c. 
The  number  density  profiles  of  the  droplei  components  (Na  +  ,  I  and 
formamide)  are  shown  in  (a)— (c),  with  line  designation  given  in  (a). 
In  (b)  and  (c)  we  show  als*o  the  difference  in  number  densities  of  the 
ions  N’sa  —  Ni  (solid  line  near  lhe  axis;  see  designation  in  (b)),  as  well 
as  the  electric  fields.  The  iniemal  field  of  lhe  droplet  (generated  by 
indueed  surface  charges)  is  shown  as  lighter  lines  and  fluctuates  around 
0.6  V/nm.  The  indueed  field  outside  the  droplet  is  given  by  squares 
(connected  by  a  line,  marked  E-field).  The  aspeci  ratios,  //,  for  lhe 
various  configurations  are  specified  in  lhe  figure 

second  or  third  interfacial  layer  (~6  A  from  the  interface)  in 
comparison  with  the  bulk.  Differences  between  Na+  and  I” 
densities  are  insignificant.  This  is  consistent  with  the  traditional 
view102  pertaining  to  the  behavior  of  the  interfacial  region  of 
salt  solutions,  and  it  contrasts  with  recent  MD  and  experimental 
findings  in  aqueous  Nal  solutions,  where  the  highest  I  densities 
were  found  in  the  outer  interfacial  layer;54*55-64  6667-102  we  defer 
further  comments  on  this  issue  to  our  discussion  in  section  5. 

In  Figure  13b  we  display  the  density  profiles  for  the  droplet 
shown  in  Figure  1  lb.  The  interfacial  charge  density  has  clearly 
increased  compared  to  the  E  ~  0  case  (Figure  1 3a).  Near  the 
positive  end  (right  in  Figure  13b),  there  is  a  slight  preference 
for  positive  ions  with  ion  densities  of  both  polarities  peaked  at 
the  interfacial  layer.  At  the  negative  droplet  end,  there  is  a  clear 
preference  for  negative  ions  in  the  interfacial  layer.  The  droplet 
electric  field  strength  (due  to  the  droplet  charges)  is  highest 
outside  the  droplet  (beyond  the  droplet  ends)  with  total  fields 
(including  the  applied  field)  exceeding  1.5  V/nm.  Interestingly, 
the  maxima  are  displaced  away  from  the  maxima  in  the  ion 
density,  suggesting  that  the  aligned  dipoles  of  the  solvent 
molecules  contribute  considerably  to  the  electric  field  in  the 
immediate  vicinity  of  the  ends  of  the  droplet.  Indeed,  when  the 
field  produced  by  the  droplets  is  decomposed  into  contributions 
from  the  ions  and  the  formamide  molecules,  we  find  that  both 
make  significant  contributions.  It  is  also  of  interest  to  note  that 
the  axial  electric  field  of  the  pure  formamide  droplet  (placed  in 
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an  external  field  of  0.5  V/nm;  see  Figure  10)  is  essentially  the 
same  as  that  shown  in  Figure  1 3b  for  the  salt  solution  droplet. 
In  the  pure  formamide  droplet  the  total  contribution  to  the 
droplet  polarization  field  comes  from  the  aligned  dipoles  of  the 
formamide  molecules.  From  the  fact  that  the  fields  in  the  exterior 
of  the  droplet  have  similar  magnitudes  in  the  pure  and  salt 
solutions  droplets  (where  in  the  latter  both  the  molecular  diploes 
and  the  dissolved  salt  ions  contribute  to  the  total  field),  we 
conclude  that  the  larger  contribution  of  the  molecular  dipoles 
in  establishing  the  field  in  the  salt- free  formamide  droplet  reflects 
a  higher  degree  of  dipole  alignment  in  such  droplets  compared 
to  that  in  the  salt  solution  ones  where  the  dissolved  ions  frustrate 
the  molecular  alignment  along  the  field  direction.  Direct  analysis 
of  the  dipole  alignments  in  pure  formamide  and  salt  solution 
droplets  verifies  the  above  conclusion.  We  also  note  that  the 
field  strengths  (of  the  order  of  1  V/nm)  originating  from  droplet 
polarization,  are  established  very  shortly  after  application  of 
the  external  field  and  before  the  occurrence  of  any  significant 
field-induced  elongation. 

The  field  produced  by  the  charges  of  the  constituents  of  the 
droplet  vanishes  asymptotically  far  from  the  droplet,  so  that  the 
total  summed  field  converges  to  the  external  constant  field  of 
0.625  V/nm  The  negative  field  inside  the  droplet  is  of  the  same 
nature  as  the  field  within  a  dielectric  slab  or  sphere  subjected 
to  an  external  field  (see  dotted  line  in  Figure  13b,c).  Charges 
of  opposite  polarity  are  induced  on  opposing  sides  of  the 
dielectric  and  this  creates  an  internal  field  that  opposes  the 
external  field. 

In  the  present  case,  the  average  internally  induced  field  is 
0.58  V/nm,  corresponding  to  a  small  total  internal  field,  E 
=  0.625—0.58  ~  0.04  V/nm.  Such  small  internal  fields  are 
typical  of  dielectric  materials  responding  to  external  electric 
fields.  The  high  sensitivity  of  computed  dielectric  properties  to 
small  variations  in  this  internal  field  necessitates  the  accumula¬ 
tion  of  sufficient  data  that  would  permit  a  meaningful  statistical 
analysis.  The  highly  oscillatory  nature  of  the  field  within  the 
drop  seen  in  our  simulations  is  due  to  limited  statistics  and  the 
fact  that  the  computed  field  reflects  a  nonuniformity  of  the  ionic 
distribution  that  does  not  self-average  over  the  averaging  interval 
of  — '3(K)  ps  used  in  these  simulations. 

4.3,3.  Charge  Emission  Statistics.  Statistics  on  charge  emis¬ 
sion  is  collected  during  the  simulation  by  recording  the  passage 
of  all  ejected  clusters  through  the  aforementioned  “detection 
planes"  placed  at  ±200  nm  along  the  z-axis  The  atomic  position 
and  momentum  data  of  each  cluster  are  also  recorded  so  that 
the  evolution  of  each  individual  cluster  can  be  examined  in  a 
separate  continued  simulation  of  the  isolated  species  (after  the 
given  cluster  is  essentially  free  from  all  the  others). 

The  surface  tension  of  the  Nal  — formamide  droplets  is 
comparable  to  that  of  the  pure  formamide  droplet  and  the 
predicted  critical  field  for  instability17  discussed  earlier  (see  also 
Figure  10)  is  similar.  However,  in  the  salt-solution  case,  the 
droplet  does  not  simply  elongate  to  another  stable  configuration 
at  the  critical  field.  Instead,  the  imbalance  between  electrical 
stresses  and  capillary  forces  is  relieved  through  a  more  favorable 
route,  namely  one  involving  formation  of  a  jet  of  charged 
clusters  (at  each  end  of  the  elongated  nanodroplet).  At  the 
smaller  investigated  fields  (<1  V/nm),  the  droplets  elongate 
more  slowly  and  emit  cluster  ions  less  frequently  w  ith  the  result 
that  capillary,  Rayleigh-type,  instabilities  lead  to  break-up  of 
the  droplet  into  large  charged  segments  that  are  accelerated  out 
of  the  computational  cell  before  the  number  of  cluster  ions  that 
are  emitted  is  adequate  for  meaningful  statistics.  Consequently, 
when  investigating  cluster-ion  emission  statistics,  we  focus  on 
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Figure  14.  Charge  abundances  of  ion-containing  clusters  emitted  from 
a  1:8  ralio  Nal— formamide  droplet  placed  in  a  uniform  exiernal  field. 
Results  are  shown  for  external  field  strengths  of  I  0  and  1 .5  V/nm 


figure  15.  Ion  composition  of  singly  charged  fomiamide/Nal  clusters 
emiued  from  a  1:8  ratio  Nal  — formamide  droplet  placed  in  a  uniform 
exiernal  field  of  1.5  V/nm  (lop)  and  1.0  V/nm  (boitom). 

simulations  of  droplets  at  the  higher  field  strengths.  The 
improved  statistics  associated  with  the  higher  charge  fluxes 
allows  for  a  better  comparison  with  the  mass  spectrometric 
measurements. 

Under  the  influence  of  an  electric  field  of  1.5  V/nm,  the 
droplet  begins  to  deform  along  the  axial  field  direction  about 
10—20  ps  prior  to  charge  emission  of  positively  and  negatively 
charged  solvated  ions  from  opposite  ends  of  the  drop.  Over  the 
course  of  about  1  ns,  the  parent  droplet  continues  to  elongate 
to  at  most  an  aspect  ratio  of  t]  =  8,  and  the  onset  of  charge 
emission  is  observed.  The  continuous  stream  of  charged  particles 
prevents  elongation  to  higher  aspect  ratios  such  as  the  one 
observed  for  the  pure  formamide  droplet  (see  Figure  10).  The 
simulations  demonstrate  that  charge  emission  happens  from 
multiple  sites  on  each  end  of  the  drop.  This  is  consistent  with 
a  stressed,  multicone  electrospray  regime  discussed  earlier.  We 
remark  here  that,  because  of  the  early  emission  of  charged 
particles  in  simulations  employing  a  high  electric  field  (1.5 
V/nm),  most  of  our  analysis  focused  on  lower  applied  fields 
which  enable  an  improved  comparison  with  the  experiments, 
particularly  on  issues  pertaining  to  the  dynamics  of  the  cluster 
ion  emission  and  droplet  formation. 

The  charge  distribution  of  cluster  ions  originating  from 
formamide  droplets  containing  29%  (by  weight  corresponding 
to  1 :8  mol  ratio)  of  dissolved  Nal  is  shown  in  Figure  14.  Results 
are  presented  here  for  two  separate  simulations  having  externally 
applied  electric  field  strengths  of  1.0  and  1.5  V/nm.  The 
abundances  of  charge  in  both  of  these  simulations  are  similar, 
showing  that  the  majority  of  the  observed  clusters  have  ±  1 
charge. 

More  detailed  information  on  the  ionic  composition  of  the 
most  frequently  observed  clusters  of  charge  ±1  is  shown  in 
Figure  15.  In  the  case  of  positive  singly  charged  clusters,  the 
fraction  of  clusters  with  ionic  composition  Naf,  Na:+1  ,  etc 
and  for  negatively  charged  clusters,  the  fraction  corresponding 
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Figure  16.  Solvation  number  distribution  of  single-ion  clusters  emitted 
from  a  1  8  ratio  Nal— formamide  droplet  under  a  uniform  external  field 
of  1.5  V/nm.  MD  results  corrected  for  30  //s  flight  time  (adjusted  MD) 
arc  also  shown  (see  text  for  details). 

to  1  ,  Na^h-,  etc,  are  shown  regardless  of  the  number  of 
accompanying  adsorbed  (solvating)  formamide  molecules. 
Though  multi-ion  clusters  are  formed,  the  primary  contribution 
is  from  single-ion  clusters  (Na+  and  1“  with  varying  numbers 
of  adsorbed  formamide  molecules).  Furthermore,  the  multi-ion 
cluster  fraction  increases  with  field  This  is  consistent  with  the 
present  experimental  single  cone-jet  observations  at  lower 
representative  field  strengths,  where  the  multi-ion  fraction  due 
to  Na2+1  and  Na+li-  solvated  ions,  is  very  small  (see  Figures 
4  and  5). 

In  Figure  16  we  compare  the  solvation  number  distribution 
of  the  single-ion  positive  and  negative  cluster  ions  compiled 
from  the  concentrated  ( 1:8  mol  ratio),  1.5  V/nm  MD  simulation 
to  those  observed  in  the  concentrated  electrospray  experiment 
Also  shown  are  MD  results  adjusted  for  a  ~30  //s  flight  time 
using  the  calculated  survival  fractions  from  Tables  3  and  4  for 
30()  K.  For  Na+  and  1  solvation  numbers  exceeding  6  and  5, 
respectively,  a  zero  survival  fraction  is  assumed,  but  for  n  <  4 
and  3,  respectively,  a  100%  survival  is  implied.  The  Na+  clusters 
have  predominantly  4-5  formamides,  and  the  adjusted  cluster 
distribution  compares  well  with  the  experimental  abundances. 
The  clusters  exhibit  a  much  broader  distribution  with  some 
preference  toward  lower  solvation  numbers.  The  adjusted  MD 
statistics  are  again  consistent  with  the  experimental  results 
(Figures  4  and  5)  for  a  single  cone-jet  electrospray  of  the  same 
concentrated  solution. 

5.  Discussion 

The  close  correspondence  that  we  found  between  the  data 
measured  in  the  vacuum  electrospray  experiments  and  the  results 
recorded  in  our  theoretical  MD  simulations  performed  on  liquid 
drops  (made  of  pure  formamide  or  Nal— formamide  solutions) 
demonstrates  that  such  coordinated,  complementary,  and  supple¬ 
mentary  studies  form  a  most  valuable  methodology,  and  an 
invaluable  source  of  information  and  deep  insights  pertaining 
to  the  complex  physics  and  operational  control  of  electrosprays, 
nanojets,  and  charged  droplets  under  the  influence  of  high 
electric  fields  Having  presented  in  the  previous  sections  our 
laboratory  measurements  and  theoretical  simulations,  we  discuss 
now  certain  issues  that  have  emerged  through  the  analysis  of 
the  results  of  our  investigations. 

The  characteristic  parameters  of  our  experiments  (see  Table 
1 )  demonstrate  that  the  cone-jet  generated  in  these  experiments 
operates  in  the  nanoscale  jet  regime  with  sufficiently  high  fields 
generated  normal  to  the  liquid  surface  to  induce  effective  field 


evaporation  of  ions.  The  observed  average  mass-to-charge  ratios 
and  the  significant  cluster-ion  currents  are  signatures  of  a  mixed 
droplet/ion  emission  regime  where  a  significant  amount  of 
current  is  produced  by  ion  held  evaporation  from  the  cone-jet 
surface,  whereas  the  mass  flow  is  dominated  by  the  droplets 
Although  pure  ion  emission  would  result  in  average  mass-to- 
charge  ratios  below  l(XX)  amu,  typical  mixed  droplet/ion 
emission  with  weak  (“background”)  ion  currents,  possibly  due 
to  droplet  ion  emission,  yields  average  mass-to-charge  ratios 
near  and  above47  10000.  The  present  average  mass-to-charge 
ratios,  which  are  below  l(XXX)  (but  well  over  1000),  point  to 
an  intermediate  regime  with  comparable  contributions  from  ions 
and  charged  droplet  currents.  This  conclusion  is  further  con¬ 
sistent  with  the  estimated  values  for  the  maximum  surface- 
normal  fields,  Emax  (between  1  4  and  1.7  V/nm),  and  for  the 
critical  radii,  /?max  (between  2  and  3  nm),  derived  from  the  axial 
emission  currents  and  mass  flows,  thus  confirming  predictions 
obtained  from  the  theory  of  cone-jet  emission.44-41  The  deduced 
jet  diameter  of  about  5  nm  is  comparable  to  that  of  the  simulated 
droplets,  supporting  the  comparison  between  the  present  elec¬ 
trospray  experiment  and  the  MD  simulations. 

Decreasing  the  salt  concentration  of  the  Nal -formamide 
solution  leads  to  a  reduction  of  the  maximum  surface-normal 
field  strength  and  an  increase  of  the  associated  radius  of 
curvature.  This  is  consistent  with  eqs  6  and  7,  The  expected 
increase  in  mlq,  however,  is  found  to  be  significant  only  for 
negatively  charged  emitted  particles.  This  somewhat  surprising 
result  can  be  explained  by  the  substantial  decrease  in  charged 
particle  beam  divergence  upon  dilution,  as  reflected  by  the  ratio 
of  the  Faraday  Cup  current  with  respect  to  the  total  emitted 
current  As  angular  distribution  measurements  on  ionic  liquid 
cone-jets  operated  in  a  nano  regime  have  shown,51 52  the  mfq 
distnbution  can  be  nonuniform  with  respect  to  the  emission  polar 
angle  Thus,  the  derived  average  mlq  value  is  only  representative 
of  charge  emission  captured  at  small  angles  with  respect  to  the 
capillary  axis  The  nonuniformity  reported  by  Chiu  et  al.  for 
the  ionic  liquid,  [Emim]Ilm|,M-52  is  consistent  with  predictions 
that  the  surface  electric  field  is  highest  at  the  cone-jct  neck, 
where  field-evaporated  ions  are  driven  away  from  the  cone-jet 
axis.41-44  The  respective  field-evaporated  ions  are,  therefore, 
preferably  emitted  at  large  angles  whereas  droplets,  emitted  from 
the  jet  tip,  have  narrow  angular  distributions  centered  on  axis. 
This  interpretation  was  further  confirmed  for  the  ionic  liquid 
by  the  measured  energy  distributions  demonstrating  that  the  ions 
peak  at  energies  nearly  equal  to  the  emitter  potential,  whereas 
the  droplets  appear  at  lower  potentials.  Consequently,  if  similar 
field  evaporation  dynamics  is  operative  in  the  present  concen¬ 
trated  Nal— formamide  solution,  a  significant  fraction  of  the  ion 
current  is  not  detected  by  the  Faraday  Cup,  and  the  average 
value  of  mlq  for  the  entire  jet  is  considerably  higher  for  cases 
with  substantial  divergence  of  the  ionic  component  in  compari¬ 
son  with  the  droplet  component.  This  interpretation,  however, 
is  contradicted  by  the  fact  that  the  present  retardation  potential 
measurements  (Figure  7)  show  that  only  a  small  fraction  of  ions 
are  produced  at  energies  exceeding  the  droplet  energy,  sug¬ 
gesting  that  most  ions  are  formed  in  a  region  of  the  neck 
downstream  from  where  ohmic  losses  have  occurred.  Similar 
retarding  potential  results  have  been  derived  by  Gamero-Castano 
for  concentrated  Nal— formamide  at  low  flow  rates.44  Conse¬ 
quently,  the  field  evaporation  dynamics  of  the  present  solution 
may  not  be  equivalent  to  that  observed  for  the  ionic  liquid 
[Emimfllm].  This  does  not  discount  the  possibility  that  ion 
emission  is  more  divergent  than  the  droplet  spray. 
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In  addition  to  field  evaporation  from  the  jet,  ions  can  he 
emitted  from  droplets  as  described  by  the  mechanism  of  lribame 
and  Thomson24  for  electrosprays  of  volatile  solvents  in  an 
ambient  atmosphere,  This  requires  that  the  surface  electric  field 
strength,  given  by  eq  2,  exceeds  the  ion  evaporation  (desorption) 
limit.  However,  this  mechanism  can  only  lead  to  a  current  that 
is  a  small  fraction  of  the  total  droplet  current  because  the 
surface-normal  field  strength  is  expected  to  decline  below  the 
field  evaporation  limit  following  the  emission  of  only  a  few 
ions,41  particularly  for  nonvolatile  droplets  of  nanometer 
dimension. 

Further  insights  were  obtained  here  through  comparative  mass 
spectrometric  studies  of  concentrated  and  dilute  Nal— formamide 
solutions.  Although  the  distribution  of  Na+  solvated  ions  is  very 
similar  at  the  two  concentrations,  apart  from  a  greater  number 
of  trace  ions  at  high  solvation  numbers,  the  l~  solvated  ion 
distribution  broadens  considerably  and  shifts  to  higher  solvation 
numbers  upon  dilution.  This  behavior  is  consistent  with  recent 
theoretical  and  experimental  investigations  of  aqueous  salt 
solutions,53  57.60  67.102  where  the  observed  preference  at  the 
interface  for  more  polarizable  anions  such  as  I  was  explained 
and  correlated,  in  part,  with  the  broader  range  of  solvation 
numbers  that  the  anion  can  accommodate,  in  comparison  with 
less  polarizable  anions  and  cations,  which  prefer  more  rigid, 
hydrogen-bonded,  solvation  shells.  Thus,  as  the  concentration 
decreases,  more  solvent  molecules  are  available  to  bond  with 
1  ,  whereas  Na+  maintains  a  tightly  bound  structured  solvation 
shell. 

TTie  fraction  of  Na+  ions  with  solvation  number  n  >  4  drops 
rapidly.  Indeed,  for  such  solvated  clusters  room  temperature 
thermal  instability  is  predicted;  the  present  quantum  chemical 
results  (see  Table  2)  predict  that  the  solvated  ions  become 
thermally  unstable  at  room  temperature  above  n  =  4.  The  1 
solvated  clusters  are  predicted  to  be  thermally  unstable  for  n  > 
3  Results  from  ab  initio  calculations  provide  a  lower  bound 
for  the  largest  value  of  n  that  provides  a  stable  cluster  at  room 
temperature.  The  calculations  assume  that  all  the  energy  that  is 
availahle  to  the  system  is  channeled  into  a  dissociative  mode. 
However,  deviations  from  the  observed  experimental  values  of 
n  suggest  that  the  energy  is  distributed  among  several  available 
modes  whose  number  increases  with  n.  Overall,  the  ab  initio 
and  RRKM  calculations  do  not  preclude  the  existence  of  I 
solvated  cluster  with  n  >  3.  In  fact,  the  larger  number  of  ions 
with  high,  thermally  unstable,  solvation  numbers  observed  upon 
dilution,  suggests  that  a  significant  number  of  cluster  ions  are 
emitted  with  large,  unstable,  solvation  numbers,  which  subse¬ 
quently  lose  (via  evaporation)  solvating  ligands,  resulting  in 
evaporative  cooling  of  the  product  cluster  ions.  This  effect  is 
even  more  pronounced  for  the  1  cluster  ions,  to  the  point  that 
the  most  abundant  solvation  numher  lies  heyond  the  present 
predictions  of  thermal  stability  at  a  liquid  temperature  of  298 
K  (Table  2). 

Similar  behavior  was  observed  in  our  experiments  on  KI 
solutions,  where  the  K+  solvated  ions  also  have  maximum 
ahundances  near  the  predicted  limit  of  thermal  stability.  The 
main  difference  between  the  KI  and  Nal  solutions  are  the 
survival  of  complex  ions  (of  both  polarities)  upon  dilution  for 
both  ion  polarities,  observed  for  the  KI  case.  The  KI  experiments 
were  mired  by  crystallization  of  the  solute  and,  therefore,  had 
to  be  conducted  at  significantly  lower  molar  concentrations.  The 
higher  propensity  for  forming  complex  ions  [X±(KI)- 
(HCONH:),,,  X  =  K^  or  I-],  despite  the  lower  concentration  in 
comparison  with  the  Nal  solution,  may  be  a  manifestation  of 
the  lower  solubility  of  the  KI  solute. 


The  snapshots  of  the  droplet  configurations  recorded  in  the 
simulation  of  a  Nal— formamide  solution  droplet  with  1:16 
(Nal— formamide)  ratio  (Figure  1 1)  under  the  influence  of  an 
electric  field  of  0.625  V/nm,  reveal  a  numher  of  features  that 
are  in  direct  correspondence  with  the  experimental  observations. 
In  particular,  in  both  cases  we  observe  mixed  droplet  formation 
and  field-induced  ion  evaporation  (desorption)  events.  Formation 
of  nanojets  is  clearly  observed  in  the  simulations  (see  Figure 
1  Id  where  a  jet  with  a  diameter  of  about  2  nm  is  shown), 
However,  we  remark  that  the  jet  emerging  in  our  simulation  is 
a  transient  state  whose  further  development  prior  to  particle 
emission  and  breakup  into  droplets  may  be  limited  by  the  finite 
size  of  the  simulated  drop.  In  the  experiment,  on  the  other  hand, 
a  continuous  (unlimited)  source  of  liquid  is  present,  thus 
allowing  the  development  of  steady-state  jet  (flow)  conditions. 
Therefore,  to  allow  a  detailed  theoretical  study  of  the  breakup 
processes,  simulations  at  lower  applied  fields  were  favorable, 
Interestingly,  in  the  simulations,  the  majority  of  cluster  ions 
are  found  to  be  emitted  from  the  tip  of  the  jet  rather  than  from 
the  neck  region  next  to  the  Taylor  cone  (for  geometrical 
definitions  see  Figures  1  and  2).  This  finding  is  consistent  with 
the  experimental  energy  distributions  of  the  solvated  ions,  which 
demonstrate  that  indeed  most  ions  are  emitted  closer  to  the  jet 
region,  that  is,  beyond  the  cone-neck  region  where  ohmic  losses 
occur.  This  observation  is  also  consistent  with  continuum 
electrohydrodynamic  predictions  of  cluster-ion  evaporation  at 
surface  regions  of  high  curvature  and  therefore  maximum 
surface  electric  field  strengths,  which  may  he  the  cone-neck 
region,  the  breakup  region  of  the  jet  (usually  near  the  tip  of  the 
jet),  or  the  emitted  charged  droplets.41*43  In  the  nanoscale  jets 
observed  in  this  study,  the  regions  of  highest  spatial  curvature 
are  at  the  ends  of  the  jets  where  nascent  drops  are  either  forming 
or  have  just  detached 

Because  of  the  small  dimensions  of  the  jets  formed  in  the 
MD  simulations,  as  well  as  deduced  from  analysis  of  the 
experimental  measurements,  one  may  indeed  expect  that  a 
continuum  approach  would  fail  to  provide  an  adequate  descrip¬ 
tion  of  the  evolution  and  breakup  processes  of  the  complex 
liquid  system  on  which  we  focus  here,  From  inspection  of  the 
droplet  configurations  displayed  in  Figure  11,  we  observe  the 
occurrence  of  spatial  and  temporal  fluctuations  of  the  ion  density 
that  are  expected  to  influence  in  a  significant  manner  the  electric 
field  distribution  and  the  consequent  stability  and  evaporation/ 
fragmetation  dynamics  of  the  system  Indeed,  we  have  noted 
that  droplet  formation  events  (location  and  timing)  correlate  with 
instances  of  ion  bunching  in  the  jet  region.  These  observations 
are  likely  to  be  of  particular  relevance  to  electrolyte  solutions 
where  bunching  of  the  solute  ions  (sodium  and  iodine  ions  in 
our  case)  occurs;  indeed,  in  the  case  of  an  ionic  liquid  (where 
ions  are  distributed  uniformly,  unlike  our  case  here  of  salt 
solutions  where  the  solute  salt  ions  are  dissolved  in  a  polar 
solvent),  certain  characteristics  that  indicate  emission  from  the 
neck  region  have  been  observed,51*52  even  though  the  jet 
diameter  in  that  system  is  estimated  to  be  similar  to  that  found 
in  the  present  study.  Additionally,  stochastic  thermal  effects  play 
a  much  increased  role  for  fluid  jets  of  nanoscale  dimensions,69 
modifying  the  kinetics  of  fluid  pinch-off  from  the  predictions 
of  a  purely  deterministic  continuum  approach.  Formulation  of 
electrohydrodynamic  theoretical  treatments  of  instabilities  and 
fragmentation,  with  the  inclusion  of  fluctuations  due  to  inho¬ 
mogeneities  in  the  solvated  ionic  charge  distribution  in  the  fluid 
as  well  as  size-dependent  thermal  noise,  present  a  formidable 
future  challenge. 


‘Feature ‘Article 


J .  Phys.  Chem .  4,  Vol.  112 ,  40,  2^  9645 


An  important  prediction  obtained  from  the  MD  simulations 
concerns  a  crossover  behavior  as  a  function  of  the  size  of 
the  emitted  species,  which  distinguishes  between  two  distict 
regimes:  one  corresponding  to  solvated  single-ion  clusters 
with  radii  R  <  1  nm,  and  the  second  one  corresponding  to 
multi-ion  droplets  that  emerge  for  larger  size  fragments,  R 
>  1  nm.  Such  crossover  is  indeed  found  experimentally,  and 
it  may  he  explained  through  electrostatic  considerations  (see 
Figure  12).  This  crossover  reflects  stronger  bonding  in  the 
smaller,  single-ion  solvated  clusters  due  to  direct  ion-dipole 
interactions  which  correlate  with  the  ability  to  sustain  larger 
values  of  charge-to-mass  ratios  compared  to  larger  size 
clusters  (see  Figure  12b). 

Emission  of  clusters  that  include  ion  pairs  was  detected 
in  both  the  experiments  and  the  simulations.  Though  only 
species  containing  a  single  ion  pair  were  observed  in  the 
experiments  for  concentrated  Nal  formamide  solutions,  and 
none  for  the  dilute  solution,  higher  fractions  of  such  ions 
with  species  containing  two  or  more  ion  pairs  were  found  in 
the  simulation  (Figure  15).  Note  that  the  respective  radii  of 
the  complex  ions  containing  ion-pairs  exceed  1  nm,  and 
exhibit  electrostatic  properties  of  droplets  according  to  the 
analysis  of  Figure  12.  We  have  not  conducted  an  RRKM 
analysis  and  thermochemical  analysis  for  these  species  to 
assess  whether  such  complex  ions  decay  faster  than  the 
experimental  time  scale.  Similarly,  no  multiply  charged  ions 
were  detected  in  the  experiments,  although  they  were  found 
in  the  simulation  (Figure  14).  Higher  field  strengths  do  not 
appear  to  be  the  source  of  these  findings  because  an  increase 
of  the  field  in  the  simulations  from  1  to  1.5  V/nm  did  not 
affect  the  observed  trends  (see  Figures  14  and  15).  We  note, 
however,  that  earlier  vacuum  electrospray  studies  conducted 
with  glycerol  solutions  observed  a  multitude  of  solvated  ions 
including  metastable  ions  and  multiply  charged  species.3,31*32103 

In  general,  the  emitted  cluster  ions  recorded  in  the  MD 
simulations  contain  more  nascent  formamide  molecules  (that 
is,  they  are  characterized  by  larger  solvation  numhers)  than  the 
clusters  observed  in  the  experiments  where  typical  flight  times 
(that  is,  the  elapsed  time  between  the  instant  of  cluster  generation 
and  the  moment  of  its  detection)  are  of  the  order  of  30  //s.  The 
recorded  occurrence  of  larger  solvated  ion  clusters  is  likely  to 
be  correlated  with  the  propensity  of  the  solvent  (formamide)  to 
form  hydrogen  bonds  (see  Figures  1 1c— e).  Because  the  largest 
solvated  ions  observed  in  the  MD  simulation  are  thermally 
unstable,  evaporative  loss  of  formamide  molecules  over  longer 
flight  times  can  effectively  lower  the  average  solvation  numbers 
to  values  close  to  those  found  experimentally  (see  Figure  16). 
The  observation  of  larger  solvated  ions  in  the  simulations 
suggests  that  such  larger  cluster  ions  are  indeed  generated  in 
the  laboratory  experiment,  and  that  subsequently  these  large, 
thermally  (room  temperature)  unstable  clusters  reduce  their  size 
via  a  cascade  of  evaporative  cooling  processes  that  occur  along 
the  ion  trajectories  to  the  detector.  Indeed,  separate  simulations 
on  large  thermally  unstable  cluster  ions  show  that  the  loss  of  a 
formamide  molecule  leads  to  an  approximate  reduction  of  about 
50  K  in  the  cluster-ion  temperature.  As  evident  from  the  data 
shown  in  Table  3,  a  50  K  change  in  temperature  has  a  rather 
significant  effect  on  the  stability  of  the  cluster  ions. 

Although  the  comparison  with  the  present  experimental 
solvated  Na+  ion  distributions,  predicted  thermal  stabilities,  and 
RRKM  decay  rates  suggest  a  surface  temperature  below  350 
K,  a  200  V  drop  between  the  capillary  potential  and  the  observed 
ion  appearance  potential  indicates  the  occurrence  of  ohmic 
heating  in  the  neck  region  of  the  cone-jet.  Thus,  a  slightly 


elevated  temperature  at  the  point  of  ion  field  evaporation  could 
be  expected.  We  also  note  here  that  in  the  MD  simulations  we 
observed  an  increase  in  the  nascent  cluster-ion  temperature  from 
310  K  to  about  400  K  from  the  beginning  to  the  end  of  the 
simulation.  This  can  be  attributed  to  the  fact  that  the  simulations 
were  conducted  at  a  relatively  high  external  field,  corresponding 
to  stressed  conditions. 

A  rather  convincing  illustration  that  the  emitted  charged 
droplets  recorded  in  the  simulations  are  consistent  with  the 
values  deduced  from  the  experimental  data  for  the  characteristic 
maximal  electric  field,  Emax,  and  the  corresponding  local 
curvature  (determined  by  the  local  radius,  /?max,  of  the  flowing 
fluid  structure),  can  be  found  in  Figure  12.  We  find  that  the 
surface  fields  of  the  droplets  are  close  to  the  experimental, 
polarity  averaged,  value  of  Em ax,  which  is  found  to  be  close  in 
value  to  the  interfacial  field  determined  from  the  simulations 
(see  Figures  12a  and  13b,c).  Additionally,  the  charge-to-mass 
ratio  associated  with  the  experimental  maximal  value  of  the 
electric  held,  £max  (eq  12),  is  found  to  agree  well  with  the  values 
derived  for  the  emitted  charged  droplets  (see  Figure  I2h). 
Finally,  the  observed  trend  in  the  emitted  droplet  charge  (as  a 
fraction  of  the  predicted  Rayleigh  limit)  is  found  to  be  consistent 
with  the  upper  limit  given  by  eq  13  when  using  the  experimental 
value  of  £max  (see  Figure  12c). 

It  is  pertinent  to  remark  here  on  the  small  asymmetry  between 
the  calculated  negative  and  positive  ion  interfacial  density 
profiles  (see  Figure  13),  particularly  in  light  of  recent  interest 
in  the  interfacial  properties  of  aqueous  solutions.53-67,102  Upon 
application  of  an  electric  held,  the  ion  density  at  the  interface 
increases  (see  Figure  13b).  At  field-free  conditions  no  discemihle 
difference  is  found  between  the  negative  and  positive  ion 
interfacial  densities.  However,  when  the  simulated  droplet  is 
placed  in  a  uniform  electric  field,  it  exhibits  a  preference  for 
negative  ions  to  be  found  near  the  negative  polarity  interface, 
whereas  positive  ions  are  preferred  at  the  positive  end  The 
derived  solution-induced  electric  field  strength  normal  to  the 
surface  at  the  negative  end  exceeds  1.0  V/nm,  corresponding 
to  a  total  field  strength  of  about  1.6  V/nm  when  including  the 
external  field.  This  value  is  again  in  excellent  agreement  with 
the  derived  values  for  the  concentrated  solution  (see  value  in 
Table  1,  obtained  with  the  use  of  eq  6),  The  field  is  slightly 
lower  at  the  positive  end.  The  large  surface  helds  at  the  polar 
ends  of  the  drop  are  not  solely  a  reflection  of  the  ionic  charge 
density  hut  are  attributable  to  both  aligned  formamide  dipoles 
as  well  as  field-induced  polarization  and  separation  of  the  ionic 
charges.  Similar  surface  fields  are  deduced  for  the  pure 
formamide  droplet  (see  Figure  10)  where,  as  noted  above,  the 
lack  of  ionic  charge  is  compensated  by  an  increased  field  due 
to  a  higher  formamide  dipole-field  alignment.  As  aforemen¬ 
tioned,  the  above  indicates  that  in  solution  droplets  the  presence 
of  solvated  ions  frustrates  the  alignment  of  the  dipolar  solvent 
molecules.  We  conclude  that  the  contribution  of  the  solvent 
(particularly  molecular  orientational  polarization)  to  the  inter¬ 
facial  field  must  be  included  in  a  proper  description  of  the 
dynamics  of  electrosprays  (emanating  from  from  jets  or  droplets) 
of  salt  solutions,  particularly  for  systems  of  nanoscale  dimen¬ 
sions  (as  well  as  in  discussion  of  the  microscopic  aspects  of 
larger  ones).  Closer  inspection  of  the  molecular  configurations 
in  the  neck  region  (see,  e.g..  Figure  I  1 )  reveals  that  the  dipoles 
of  the  interfacial  solvent  molecules  near  the  neck  are  aligned 
by  the  field  more  in  the  direction  of  the  flow  than  normal  to  the 
surface.  Consequently,  the  field  evaporation  force  normal  to  the 
surface  is  reduced,  thereby  providing  an  explanation  for  the  lack 
of  ion  emission  at  the  neck. 
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The  observed  preference  of  negative  ions  at  the  interface, 
under  the  influence  of  an  external  electric  field,  ts  not  as 
pronounced  as  the  corresponding  findings  at  aqueous  interfaces 
under  field-free  conditions,  where  large  polarizable  anions  were 
particularly  preferred  over  cations  in  the  interface  layer.54-55 
From  the  fact  that  the  simulated  field-free  Nal— formamide 
droplet  exhibits  no  difference  between  anionic  and  cationic 
densities  in  the  interface  layer  (where  ion  depletion  occurs  as 
expected  from  a  solvation  energy  consideration),  we  conclude 
that  the  minor  preference  for  the  anions  to  reside  in  the 
interfacial  layer,  found  in  simulations  under  the  effect  of  an 
external  field,  is  related  to  the  interaction  of  the  ions  in  the 
cluster  with  the  applied  field.  Because  the  solvation  number 
distributions  are  broader  for  the  larger  anions,  and  because  the 
ligand  binding  energies  are  weaker  (see  Table  2),  it  may  be 
expected  that  the  propensity  of  the  formamide  molecules 
coordinated  with  the  I-  to  align  with  the  field  is  higher  than 
that  of  molecules  coordinated  with  the  Na+  cations,  resulting 
in  a  stronger  attraction  toward  the  high-field  region  at  the 
interface. 

Additionally,  the  lack  of  charge  preference  in  the  field- 
free  Nal  — formamide  droplet  could  also  be  attributed  to  the 
nanometer-scale  droplet  confinement  effects,  surface  curva¬ 
ture,  and  characteristics  of  the  liquid  solvent  (that  is, 
formamide  in  the  present  simulations  versus  water  in  previous 
interfacial  studies).  The  chemical  identity  factor  could  be  of 
principal  significance,  particularly  in  light  of  the  fact  that 
the  behavior  of  aqueous  solutions  is  often  discussed  in  terms 
of  hydrogen  bonding  between  water  molecules  in  clusters  of 
large  anions,  versus  symmetric  solvation  in  the  case  of  the 
cationic  systems.102  Accordingly,  it  has  been  shown  through 
MD  simulations  of  aqueous  Naf  and  Cl  clusters,  that  Na+ 
is  solvated  symmetrically  in  the  center  of  the  clusters,  and 
Cl  resides  on  the  surface  of  a  water  cluster. 104  This  was 
confirmed  experimentally  for  the  larger  halide  ions  in  water 
dimer  clusters  by  Johnson  and  co-workers. 105  The  present 
ab  initio  optimized  structures  of  formannde-solvated  ions, 
however,  show  no  evidence  for  hydrogen  bonding  with  both 
the  formamide-solvated  Na+  and  1  ions  located  close  to  the 
centers  of  the  corresponding  clusters.  This  comparison  calls 
for  similar  simulations  of  electrospray  dynamics  in  solutions 
with  volatile  solvents  such  as  those  used  in  mass  spectro- 
metric  electrospray  experiments  involving  an  electrospray 
source  operating  at  atmospheric  pressures.  The  general  jet 
dynamics  described  in  the  present  work  are  applicable  to 
volatile  solvents  in  an  atmospheric  source.  The  present  work, 
however,  demonstrates  that  irrespective  of  solvent  volatility 
and  gas-phase  pressure,  the  detailed  field-evaporation  dynam¬ 
ics  need  to  be  examined  on  a  case-by-case  basis  because  they 
arc  governed  by  complex  solvation  dynamics  at  the  liquid— air 
interface  and  its  response  to  electric  fields  with  strengths  on 
the  order  of  1  V/nm.  We  remark  here  that  to  perform  MD 
simulations  on  aqueous  drops  would  require  the  inclusion 
of  a  gaseous  atmosphere,  which  constitutes  an  added 
complexity  in  the  simulation.  Additionally,  comparisons  with 
experimental  data  on  nascent  ion  distributions,  enabled  here 
(in  part)  by  the  low-vapor  pressure  of  the  solvent  (forma¬ 
mide),  are  not  available  for  aqueous  solutions. 

6.  Summary 

In  this  paper  we  presented  laboratory  experiments  and 
theoretical  investigations  aiming  at  characterizing  and  under¬ 
standing  the  atomic-scale  structure,  energetics,  stability,  dynam¬ 
ics,  and  fragmentation  processes  of  complex  liquids  under  the 


influence  of  high  electric  fields  The  present  MD  simulations 
of  drops  of  Nal— formamide  solutions  placed  in  a  uniform  field 
reproduce  faithfully  the  main  salient  experimental  features 
observed  in  vacuum  electrospray  ionization  experiments  of 
solutions  of  similar  concentration,  operated  in  a  single  cone-jet 
mode  with  jet  dimensions  comparable  to  the  drop  size.  At  these 
conditions,  Taylor  cones  spray  charged  particles  in  a  mixed 
solvated  ion/droplet  regime,  The  MD  simulation  of  a  droplet 
subjected  to  a  uniform  electric  field  of  0.625  V/nm  shows  the 
occurrence  of  droplet  fission  and  field  ion  evaporation  processes. 
The  simulation  predicts  that  ion  emission  occurs  at  the  breakup 
region  of  jets  (sec  Figures  2  and  1  ld,e)  where  nascent  charged 
droplets  form  and  get  emitted,  as  well  as  at  the  tips  of  elongated 
charged  droplets  and  the  nanojets  that  emanate  from  them.  This 
is  consistent  with  studies  that  correlate  enhanced  cluster-ion 
emission  with  the  regions  of  highest  spatial  curvature.41-43  The 
results  from  the  simulation  are  consistent  with  the  present 
experimental  findings  that  most  of  the  field  evaporated  ions  are 
formed  with  energies  near,  or  lower,  than  those  of  the  charged 
droplets.  The  surface  normal  field  strength  determined  from  the 
simulation  near  the  point  of  charge  emission  (for  the  concen¬ 
trated  solution  droplet)  is  determined  to  be  ^1.6  V/nm,  in 
excellent  agreement  with  the  field  derived  from  the  experimental 
observations  using  continuum  electrostatic  theory.  The  size 
dependence  of  the  surface  electric  field  and  charge- to-niass  ratios 
of  droplets  emitted  in  the  MD  simulation  are  fully  consistent 
with  the  experimentally  derived  maximum  surface  field  strength 
and  associated  curvature.  A  transition  in  the  electrostatic 
properties  of  the  emitted  clusters  as  a  function  of  size  is  observed 
at  radii  of  ~~  1  nm,  corresponding  to  a  transition  between  solvated 
single  ions  (for  R  <  Inin)  and  multi-ion  droplets  (for  R  >  lnm). 
This  crossover  is  attributed  to  the  stronger  intermolecular  forces 
and  the  consequently  higher  sustainable  charge  densities  in  the 
solvated  single  ions. 

The  simulated  and  measured  solvation  number  distributions 
for  Na+(HCONH2)„  and  I  (HCONH2),,  are  in  good  agreement. 
Though  the  Naf  cluster-ion  distributions  are  governed  primarily 
by  thermal  stability  and  are  independent  of  concentration,  the 
1  cluster-ion  distributions  broaden  with  dilution.  Similar 
behavior  is  observed  experimentally  with  KI— formamide  solu¬ 
tions  Larger  cluster  ions  that  are  thermally  unstable  at  3()0  K 
are  observed  and  attributed  to  evaporatively  cooled  species. 
Complex  ions  including  ion  pairs  are  observed  both  in  experi¬ 
ments  and  in  the  simulations.  Although  multiply  charged  ions 
are  not  observed  experimentally,  minor  quantities  of  multiply 
charged  ions  survive  the  ~15  ps  flight  time  of  the  simulation 
space. 

The  present  MD  simulations  reveal  that  fora  Nal -formamide 
droplet  at  field-free  conditions,  ion  densities  are  depleted  in 
the  interfacial  region  and  that  there  is  no  preference  of  anions 
over  cations  (see  Figure  13a).  Application  of  a  field  polarizes 
the  droplet  through  molecular  dipole  alignment  and  ion 
charge  redistribution  resulting  in  excess  positive  charge  (and 
cation  concentration)  and  excess  negative  charge  (and  anion 
concentration)  at  the  opposing  ends  of  the  droplet.  The 
resulting  electric  field  strength  at  the  exterior  surface  of  the 
droplet  is  due  to  both  the  molecular  dipole  alignment  and 
polarized  ion  distribution  and  is  of  a  magnitude  1  V/nm 
(leading  to  a  total  electric  field  of  M.5  V/nm  when  the 
externally  applied  field  is  added). 

In  summary,  the  combined  theoretical  and  experimental 
investigations  presented  in  this  paper  provide  deep  insights 
pertaining  to  the  microscopic  origins  of  the  properties  of  charged 
dielectric  fluids  (that  is,  salt  solutions  of  formamide)  with  (and 
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without)  the  influence  of  externally  applied  electric  fields.  In 
particular,  these  studies  have  unveiled  certain  molecular-scale 
aspects  of  the  instabilities  and  mechanisms  leading  to  breakup 
and  fragmentation  of  liquid  structures.  These  studies  allow  us, 
in  addition,  to  assess  the  validity  of  continuum  hydrodynamic 
and  electrohydrodynamic  theories  for  the  description  of  charged 
fluids  of  nanosale  dimensions,  and  they  may  serve  as  a  guide 
for  the  testing  and  development  of  improved  treatments,  such 
as  the  stochastic  hydrodynamic  formulation  formulated  and 
applied  recently  in  studies  of  liquid  neutral  nanojets  and 
nanobridges.69*70 
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TABLE  Al:  B3LYP/6-31G(d)  Predicted  Geometries  of 
Na(HCONH2)„%  n  =■  1-4° 


Na(HCONH2)/ 


parameter 

n  =  1 

n  =  2 

n  =  3 

n  =  4 

R(  Ol-Cl) 

1.248 

1.244 

1.232 

1.228 

R(C)  —HI) 

1.097 

1.098 

1  101 

1.103 

tf(CI-Nl) 

1.331 

1.334 

1.34 

1.344 

R(  N1-H2) 

1.012 

1.012 

1.011 

1.01 

R(  N1-H3) 

1.014 

1.014 

1.013 

1.012 

R(  02— C2) 

1.244 

1  232 

1.228 

tf(C2-H4) 

1.098 

1.101 

1.103 

R(  C2— N2) 

1 .334 

1.34 

1.344 

R(  N2— H5) 

1.012 

1.011 

1.01 

R(  N2-H6) 

1.014 

1.013 

1.012 

R(  03-C3) 

1.232 

1.228 

/f(C3-H7) 

1.101 

1.103 

R(  C3-N3) 

1.34 

1.344 

R(  N3— H8) 

1.01  1 

1.012 

R(  N3-H9) 

1.013 

1.01 

R(  04-04) 

1.228 

tf(C4-HI0) 

1.103 

R(C  4-N4) 

1 .344 

R{  N4-H11) 

1.012 

/?(N4— H  1 2) 

1.01 

a(Ol-Cl-Hl) 

1213 

1 2 1 .5 

121.8 

122 

afOl-Cl-Nl) 

124 

124.1 

124.5 

124.6 

a(Hl-Cl-Nl) 

114  6 

1  14  4 

113.7 

113.4 

a(Cl  -N 1 -H2) 

121.2 

121.3 

121.4 

121.5 

a(Cl— N1  — H3) 

120.3 

120.1 

120.1 

1  19.9 

a(H2— N 1  — H3) 

118.5 

1 18.6 

1 18.5 

1  18.7 

a(02— C2— H4) 

121.5 

121.8 

122.1 

a(02— C2— N2) 

124  1 

124.5 

124.6 

a(H4-C2-N2) 

1  14.4 

1 13.7 

1  13.4 

0UC2-N2-H5) 

121.3 

121.4 

121.5 

a(C2— N2-H6) 

120.1 

120.1 

1 19.9 

a(H5-N2-H6) 

118.6 

1 18.5 

1  18.6 

a(03-C3-H7) 

121.8 

122.1 

a(03-C3-N3) 

124.5 

124.5 

a(H7-C3-N3) 

1 13.7 

113.4 

a(C3— N3— H8) 

121.4 

119.8 

a(C3-N3-H9) 

120.1 

121.5 

a(H8-N3-H9) 

1 18.5 

118.7 

a(G4— C4— H 10) 

122.1 

a(04— C4— N4) 

124.5 

a(H  10— C4— N4) 

1  13.4 

a(C4— N4— H 1 1 ) 

1 19  8 

a(C4-N4— H 1 2) 

121.5 

a(H  1 1  — N4— 25) 

118.7 

“  Bond  lengths  in 

angstroms; 

bond  angles 

in  degrees. 

Refer 

Figure  8a— d  for  the  atomic  numbering  scheme. 


TABLE  A2:  MP2/6-311+G(d,p)  Predicted  Geometries  of 
I(HCONH2)fl  ,  n  =  1-5* 


l(HCONFh),, 

parameter  n  =  1  n  =  2  n  =  3  n  =  4  /;  =  5 


R{  Ol-Cl) 

1.23 

1.228 

1.227 

1.225 

1.222 

R(  Nl-Cl) 

1 .349 

1.35 

1.353 

1.355 

1.362 

R(  N1-H2) 

1.028 

1.022 

1.02 

1.018 

1.018 

R(  N1-H3) 

1.01 

1.011 

1.01 

1.01 

1.011 

R(C  1  — HI ) 

1.103 

1.104 

1.103 

1.104 

1.104 

R(  D2-C2) 

1.227 

1.227 

1.225 

1.224 

R(  N2— C2) 

1.355 

1.353 

1.354 

1.356 

R(  N2— H5) 

1.025 

1.02 

1.017 

1.015 

R(  N2— H6) 

1.011 

1.01 

1.011 

1.011 

R(C  2-H4) 

1.103 

1  104 

1.105 

1  105 

R(  G3-C3) 

1.227 

1.225 

1.223 

R(  N3-C3) 

1 .353 

1.355 

1.354 

/?(N3— H8) 

1.02 

1.018 

1.011 

R(  N3-H9) 

1.01 

1.01 

1.012 

/?(C3-H7) 

1  103 

1.104 

1.106 

R(  G4-C4) 

1.224 

1  221 

R(  N4-C4) 

1.361 

1  368 

R(  N4-H1  1) 

1.02 

1.021 

R(  N4-H12) 

1.01 1 

1.012 

R(C  4-H10) 

1.104 

1.103 

R(  05 -C5) 

1.224 

R(  N5-C5) 

1.356 

K(N5— HI4) 

1.017 

R(  N5-H15) 

1.01 

R(C5- HI  3) 

1.104 

a(Ol-Cl-NI) 

126 

126 

125.7 

125.5 

125.4 

a(Ol-Cl-HI) 

121.9 

121.9 

122.2 

122.3 

122.5 

a(Cl  -N 1  -H2) 

120.1 

1 20. 1 

120.1 

1 19.8 

1  1 8.5 

a(Cl  — N 1  —  H3) 

117.3 

118.5 

1 17.6 

1  17.4 

116.9 

a(N  1  — Cl— H 1 ) 

112.1 

112.1 

1  12.1 

112.2 

112.1 

a(H2  — N 1  —  H3 ) 

118.2 

117.7 

118.5 

1 18.3 

1  17.9 

a(02-C2-N2) 

125.7 

125.7 

125.8 

125.6 

a(Q2— C2-H4) 

122.2 

122  2 

122.1 

122.1 

a(C2— N2-I15) 

1 19.1 

120 

1 19.9 

1 19.8 

a(C2-N2-H6) 

1 17 

1 17.6 

1 18.6 

1  18.5 

a(N2-C2-H4) 

1 12.1 

112.1 

112.1 

1  12  2 

a(H5-N2-H6) 

1 17.9 

118  5 

117.5 

1  17.3 

a(03-C3— N3) 

125.7 

125.6 

125.9 

a(C)3— C3-H7) 

122.2 

122.2 

122 

a(C3-N3-H8) 

120 

1  19  9 

120.2 

a(C3— N3— H9) 

1 17.5 

117.6 

1 20.3 

a(N3— C3— H7) 

112.1 

112.2 

112.1 

a(H8-N3-H9) 

1 18.4 

118.4 

1  17 

a(04— C4— N4) 

125.4 

125.1 

a(O4-C4-H10) 

122.4 

122  6 

a(C4-N4-H  1 1) 

118.8 

1  17.8 

a(C4— N4-HI2) 

116.9 

1  15.8 

a(N4— C4— H  10) 

112.2 

1  12.2 

a(HI  1-N4-HI2) 

1  17.8 

117 

a(05-C5-N5) 

125.5 

a(Q5-C5-H13) 

122.3 

a(C5-N5-HI4) 

1  19  6 

a(C5-N5-H15) 

1  17  6 

a(N5— C5— HI3) 

1  12.1 

a(H14— N5  — H15) 

1  18.5 

a  Bond  lengths  in 

angstroms;  bond  angles  in 

degrees. 

Refer  to 

Figure  9a— e  for  lhe  alomic  n 

lumbering  scheme 
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Appendix 

In  this  Appendix  (Tables  Al  and  A2)  we  give  the  geometrical 
parameters  of  the  theoretically  optimized  structures,  calculated 
for  clusters  of  salt  ions  solvated  by  formamide  molecules. 
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